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ABSTRACT 



^ , The Guide Star Catalog II (GSC-II) is an all-sky database of objects derived from the un- 

$-H ■ compressed Digitized Sky Surveys that the Space Telescope Science Institute has created from 

c/3 I the Palomar and UK Schmidt survey plates and made available to the community. 

^, ' Like its predecessor (GSC-I), the GSC-II was primarily created to provide guide star infor- 

mation and observation planning support for Hubble Space Telescope. This version, however, 
is already employed at some of the ground-based new-technology telescopes such as GEMINI, 
K^ I VLT, and TNG, and will also be used to provide support for the James Webb Space Telescope 

(JWST) and Gaia space missions as well as the Large Sky Area Multi-Object Fiber Spectroscopic 
Telescope, one of the major ongoing scientific projects in China. 



^v^ ■ Two catalogs have already been extracted from the GSG-II database and released to the 



astronomical community. A magnitude-limited {Rp — 18.0) version, GSG2.2, was distributed 
(^ ■ soon after its production in 2001, while the GSG2.3 release has been available for general access 

OO ! since 2007. 

^^ ■ The GSG2.3 catalog described in this paper contains astrometry, photometry, and classifica- 

tion for 945,592,683 objects down to the magnitude limit of the plates. Positions are tied to the 
. , ■ International Celestial Reference System; for stellar sources, the all-sky average absolute error 

rS I per coordinate ranges from 0."2 to 0."28 depending on magnitude. When dealing with extended 

jrt ■ objects, astrometric errors are 20% worse in the case of galaxies and approximately a factor of 

2 worse for blended images. 

Stellar photometry is determined to 0.13-0.22 mag as a function of magnitude and photo- 
graphic passbands {Rf, Bj, In). Outside of the galactic plane, stellar classification is reliable to 
at least 90% confidence for magnitudes brighter than Rp — 19.5, and the catalog is complete to 
i?F=20. 

Subject headings: astrometry - astronomical data bases: miscellaneous - catalogs - surveys - techniques: 
image processing - photometric 



^ Space Telescope Science Institute, 3700 San Martin ence Institute, 3700 San Martin Drive, Baltimore, MD 

Drive, Baltimore, MD 21218, USA 21218, USA 

■^INAF— Osservatorio Astronomico di Torino, Strada Os- ^Steward Observatory, 9425 N. Weather Hill Dr., Tuc- 

servatorio 20, 10025 Pino Torinese, TO, Italy son, AZ 5743-5484, USA 

^Computer Sciences Corporation, Space Telescope Sci- ^Gemini Observatory Southern Operations Center, c/o 



1. Introduction 

The Guide Star Catalog II (GSC-II) is an as- 
tronomical database constructed from the scanned 
images of 9541 Palomar and UK Schmidt photo- 
graphic sky survey plates digitized at Space Tele- 
scope Science Institute (STScI). These same plate 
images are also known as the Digitized Sky Sur- 
vey (DSS), and are accessible directly from the 
STScI eb site as well as from a number of ma- 
jor astronomical data centers around the world. 
A subset of the original images, based on the 
second-epoch Palomar surveys only (also known 
as DPOSS), is availabl e separately from Caltech 
( Diorgovski et alj 120031 ). Whilst all the images 
represent over 8 terabytes of data that are archived 
at STScI, t he use of the H-transform compression 
techn ique (jWhite et al.lll992l : IWhite and Percivall 
19941 ) has enabled a reduction in data volume to 
?» 1 terabyte for distribution to data centers. 

Each individual uncompressed image was pro- 
cessed to detect the objects, and calibrations were 
obtained for each plate by polynomial modeling 
against photometric and astrometric reference cat- 
alogs. The average accuracy in position for stellar 
objects of intermediate magnitude {V < 18.5) is 
~ 0.3", whilst the corresponding astrometric pre- 
cision, on a scale of approximately 0.5 degrees, is 
close to 0.2". Photometric calibrations, carried 
out in the natural plate passband, show typical 
errors of 0.1-0.2 mag and systematic offsets lower 
than 0.1 mag. Stars are correctly classified as 
such with at least 90% confidence for magnitudes 
brighter than Rp = 19.5. 

The plate catalogs were loaded into a database 
system known as Catalog of Objects and Measured 
Para meters from All Sk y Surveys (COMPASS), 
(see iLasker et al.l Il998l ). This database was 
built using "Objectivity", a commercial Object- 
Oriented Database Management System; it consti- 
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tutes a repository for « 5 billion measurements, 
and is « 2.5 terabytes in size. Once the indi- 
vidual observations have been cross-matched to 
link different measurements of unique astronom- 
ical objects, a catalog can be exported from the 
database as a collection of FITS binary tables. 
Then, the export catalog is made available to the 
community via Web services and also provided to 
a number of data centers. Over the last few years 
there have been a number of interim releases from 
the GSC-II, primarily to consortium members for 
telescope operations. The latest public release, 
the GSC 2.3.2, contains 945,592,683 unique ob- 
jects and is ~ 170 GB in size. Details about the 
release content and access methods are available 
from the ST Scl and Osservatorio Astrononiico di 
TOrino (OATo) Web sites. 

In this paper we describe the construction, cali- 
bration, and overall quality of GSC 2.3.2 , the cur- 
rent catalog release (for simplicity, GSC 2.3 in the 
rest of this paper), which was investigated through 
internal tests and a number of external compar- 
isons to suitable literature data. 

1.1. Background 

In the age of new-technology telescopes and 
space-based missions, it is useful to recall that 
an important part of observational astronomy 
has historically been the creation of catalogs 
containing the reference and (often) target ob- 
jects required to support observing programs. 
From a historical viewpoint it is fascinating that 
astrometry — the oldest branch of astronomy — is 
so vital to the success of modern high-tech ob- 
servatories. As the technological complexity (and 
cost) of building and operating telescopes has in- 
creased enormously over the last 20 years, so has 
the effort to provide the best scientific return for 
these investments. One of the important issues 
is the optimization of observing efficiency, which 
depends on the use of proper pointing and in- 
put catalogs as well as on the access to digitized 
versions of the optical sky surveys. 

The need for a deeper, all-sky catalog was high- 
lighted early by the creation of the first Guide Star 
Catalog (GSC-I) to satisfy the pointing require- 
ments of the Hubble Space Telescope (HST). An 
in-depth description of the GSC-I catalo g may 



be found in a set of thre e papers (Lasker et al. 
19901 : iRussell et all Il990t I Jenkner et all Il990[ ). 



The GSC-I was used for HST observation plan- 
ning as well as target acquisition and tracking 
and has proved to be very reliable for its intended 
purpose. In addition, since its publication on CD- 
ROM, this catalog had become widely used by 
many ground-based telescopes to speed up the 
process of finding guide stars. 

Although the GSC-I has been used with great 
success operationally and scientifically, it became 
clear (even during its construction) that it was 
possible to improve its usefulness by addressing 
the known systematic calibration errors. Indeed, it 
was clear that an increase in scope to include mul- 
ticolor, multi-epoch data would eventually lead 
to the requirements for a new much improved 
archive, the GSC-II. In addition to the risk of dam- 
aging the UV-sensitive MAMA detectors in the 
second generation of HST instruments by O stars 
as faint as 19th magnitude, it was generally ac- 
knowledged that the proper motions of the guide 
stars during the proposed lifetime of HST would 
pose a potential problem for the accurate pointing 
of the telescope during the latter years of its op- 
eration. It was also realized that a fainter archive 
with multicolors such as the GSC-II would address 
a number of other astronomical needs such as: 
supporting adaptive optics on the next generation 
of large telescopes, remote or queue scheduling ca- 
pabilities, and improved detector instruments re- 
quiring color information of the target objects. 

In order to be prepared for this, STScI nego- 
tiated access to the original plates of the Palo- 
mar Observatory Sky Survey (POSS-II) and, in 
partnership with the Anglo-Australian Observa- 
tory (AAO), undertook the Second Epoch South- 
ern (SES) red survey with the UK Schmidt Tele- 
scope Unit (UKSTU). These new surveys, when 
combined with the plate material used in the GSC- 
I and the earlier POSS-I and UKSTU surveys, pro- 
vide the material to generate a GSC-II based on, 
at least, two epochs and three passbands. 

1.2. The GSC-II Project 

STScI and OATo first began a collaboration 
to start developing GSC-II in 1989. STScI was 
primarily motivated by telescope operations, and 
OATo was interested in the scientific applications 
of this catalog for galactic structure. 

Once the project was underway, additional re- 



sources were contributed by ESO and GEMINI 
(who wished to use GSC-II for VLT and GEM- 
INI telescope operations respectively), ST-ECF 
(as part of the NASA-ESA HST funding agree- 
ment), and the Astrophysics Division of ESA (for 
science projects). A significant fraction (~ 30%) 
of the plate processing was performed at ST-ECF 
using the same pipeline software system described 
later. 

2. Astronomical data 
2.1. Photographic Plates 

All of the survey plates scanned in this project 
were originally taken by the Palomar Schmidt tele- 
scope in California, USA, and the UK Schmidt 
telescope at Coonabarabran, Australia. Table [1] 
summarizes the main characteristics of the vari- 
ous plate material and a complete descri ption of 



the ph otographic surveys can be found in [Morgan 
(|l995l ) and references therein. 



In general, there is a strong similarity between 
the northern and southern surveys since they used 
the same photographic emulsions; however, there 
are small differences in some of the filters used, 
which result in slight differences in the transmis- 
sion curves as shown in Figure [TJ 

2.2. Astro metric Reference Catalogs 

Following the recommendation of the lAU 
XXIII General Assembly (resolution B2), the In- 
ternational Celestial Reference System (ICRS)was 
adopted to be the reference system for the GSC- 
II. At optical wavelengths, the ICRS is defined by 
the Hipparcos catalog; however, the bright Hippar- 
cos stars are heavily saturated on Schmidt plates. 
On the oth er hand, the fainter stars {V > 10) of 
th e Tycho dPerrvinan and ESAlll997l ) and Tycho- 
2 ( Hog et al.ll2000 n catalogs, do have measurable 
images on the Schmidt plates, so that each survey 
plate is directly tied to the ICRS through Tycho's 
faint end. Because the Tycho-2 catalog was not 
available until after early 2000, a subset of the 
GSC 2.3 astror netric calibrations are based on the 
ACT catalog (jUrban et al.lll998l) . Slight devia- 
tions of the ACT and Tycho-2 catalogs from the 
ICRS system are well below our measuring error. 
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Fig. 1. — Transmission curves of the photographic passbands Bj,Rp^ and I^ for the Palomar (sohd hnes) 
and AAO (dashed hne) Schmidt surveys, compared to the Johnson-Kron-Cousins BVRIc filters. 



2.3. Photometric Reference Catalogs 

The photometric cahbrators used for the GSC- 
II are the GSPC-II, GSPC-I and Tycho catalogs, 
which sample respectively the faint, intermediate 
and bright range of the photographic magnitudes. 

The GSPC-I (G uide Star Photometric Catalog 
Lasker et al.lll988l ). which was built to support the 
construction of GSC-I, contains 1477 photoelectric 
sequences with B and V standard photometry of 
9-15th mag stars at the 0.1-0.15 mag level of ac- 
curacy. 

Because of the nonlinear response of the pho- 
tographic plate, the need for fainter calibrators 
was a requirement for the construction of GSC- 
II. A lack of such calibrators in the astronomi- 
cal literature motivated the start of a decade-long 
series of observations, only now coming to com- 
ple tion, which led to the construction of GSPC - 
II (iBucciarelh et al.ll200ll : iBucciareUi et al.ll2006h . 
This is an all-sky catalog of CCD photometric se- 
quences centered in each survey plate field. Each 
sequence provides B, V and R standard photom- 
etry at the 0.05-0.1 mag precision level down to 
approximately V=19.5. The number of stars per 
sequence varies from a dozen to several hundred 



depending on the field of view of the CCD and 
the galactic latitude. At the time of GSC 2.3 cal- 
ibrations, 0.5% of the plates still did not have an 
extended photometric sequence. In these cases, 
the photometric calibration relied on existing stan- 
dards and deriving artificial points, as explained 
in Section 3.3.2. Figure [H shows the magnitude 
limit distribution of all GSPC-II calibrati on se- 
quen ces. The Tycho catalog (Perrvman and ESA 



19971 ) was used in the GSC-II photometric calibra- 
tion to stabilize the bright end of the density-to- 
intensity curve. The calibration limits of GSC 2.3 
are given in Fig. [3] which shows the distribution 
of the faintest GSCP-II reference stars available 
in the Hierarchical Triangulated Mesh (HTM) re- 
gions and used for the photometric calibration of 
the blue Bj and red Rp plates. 

3. Plate processing pipeline 
3.1. Plate digitization 

The processing of photographic plates begins 
with the use of microdcnsitometers, machines ca- 
pable of measuring photographic transmission T. 
This is the fraction of light measured as a cur- 
rent, transmitted by the portion of the developed 



emulsion as illuminated. The PDS machines used 
for the digitization of GSC-II plates measure the 
semispecular density which is defined as the light 
transmitted and diffused within an angle of < 24° . 
This is then converted into photographic density 
p through the relation 

p == - logio T. 

As T can assume all values from zero (complete 
opaque emulsion) to 1 (perfectly transparent me- 
dia), p can vary from zero to infinity; in practice, p 
usually lies in the range — 5.0. The Perkin-Elmer 
2020G PDSs originally selected for digitizing the 
STScI plate collection were first modified to cope 
with the accuracy requirements impo sed by the 
GSC-I project (see iLasker et al.lll990l) . The in- 
crease in the number of plates (5 x ) to be digitized 
for the GSC-II and DSS projects also required a 
significant increase in scanning throughput. For 
this reason, the two scanning machines available 
at STScI were refurbished with entirely new criti- 
cal subsystems to increase their speed without af- 
fecting their spatial and photometric capability; 
the improved twin digitizing devices were name d 
GAMMA 1 and GAMMA 2 (JLaidler et al.lll994l ). 

The Guide Star Automatic Measuring MA- 
chine (GAMMA) is a laser-illuminated multichan- 
nel scanning microdensitometer, modularly built 
upon the substrate of the modified PDS used in 
earlier STScI work. 

A rebuild of the x- and y-servos is based on an 
HP 5507 laser transducer system with a custom 
card to implement the PDS "1/iV" functionality. 
The light source is a spatially filtered 2 mW HeNe 
laser beam, expanded to 1 mm diameter (1/e^ of 
total energy), moved with a TeO(2) acousto-optic 
defiector (AOD), and finally imaged on the plate 
as a 42 /im Gaussian spot. This size gives low 
aliasing with the adopted 15 pvn sampling. The ta- 
ble moves in the x-direction at 100 mm s~^, while 
the AOD steps the laser beam through a small 
number of channels (typically 5) in the y-direction. 
The AOD is controlled by a frequency synthesizer 
operated from a digital signal processor (DSP) 
which cycles through a channel-to-frequency ta- 
ble, synchronized by pulses from the servo 1/iV 
logic. The DSP also removes y-position errors due 
to residual servo auto-lock disequilibrium by ap- 
plying a small frequency correction to the AOD, 
based on a linear function of the error signal from 



the y-servo. 

The light-collecting optics consists of a stan- 
dard microscope, an exit slit, an integrating sphere 
to make the response insensitive to the channel 
number, and a photomultiplier. Pixel integration 
times are about 25-50 /isec, and the data conver- 
sion is done with a 15-bit fioating point analog-to- 
digital-converter (ADC), which supplies the pho- 
tometric parameter Dij — a + bpij linearly related 
to the optical density measured at pixel i,j. 

A separate photodiode to monitor the laser 
flicker is also provided for photometric normal- 
ization. A dedicated VAXStation with IEEE-488 
and CAMAC interfaces controls each GAMMA 
machine. Separate programs operate each sub- 
system (servo, data collection, laser monitoring, 
AOD, console) and communicate with each other 
via shared memory and common event flags. This 
matching of the software and hardware architec- 
tures results in easily maintainable code. 

Nominal performance of the scanning machines 
was verified on a regular basis using several nu- 
merical methods. A sample of three areas from 
each scan was tested using a scan-line correlation 
and constrained minimization to test for scan-line 
"shear" as a function of channel. Scans where the 
shear exceeded 0.1 pixel were rejected and further 
tests and mechanical maintenance were then per- 
formed. In general the intrapass channels showed 
a very high degree of correlation and shear rarely 
exceeded 0.05 pixel. Periodic engineering scans 
were also performed on a routine basis of Max 
Levy 2 slanted and vertical rulings in order to tune 
the AOD system for optimal performance. 

3.2. Image processing 

After digitization of the Schmidt plate material 
on the GAMMA machines, each plate scan was 
processed through an image processing pipeline, 
as detailed in the following sections. 

3.2.1. Sky Background Determination 

A single value of the sky background standard 
deviation is calculated for the whole plate by sub- 
dividing it into a regular grid of 0.5 mm cells for 
the estimation of local density histograms and pix- 
els statistics. Then, the median values of each cell 
are computed, and a sequence of median filters is 
applied at different scales (8,4,2,1 cells) to produce 
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Fig. 3. — GSC 2.3 calibrated magnitude limit. Left panel: distribution of the faintest GSPC2 reference star 
available in the HTM regions and which result from the photometric calibration of the blue Bj and red 
Rf plates. Right panel: cumulative distribution of the HTM regions has a function of the faintest GPSC2 
reference stars. Note that ~90% of the sky has been calibrated with photometric sequences down to Rp > 18 
and Bj > 19.5, while for ^50% of the sky the photometric sequences attain Rp > 19.5 and Bj > 20.5. 
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Fig. 2. — Histogram of GSPC-II sequences mag 
nitude limit. 



an average sky image. The a calculation for each 
cell is found by taking the slope of the cumulative 
histogram of pixel values in the cell at the 50% 
level, to minimize the disturbance of stellar contri- 
bution. Pixels at this level are nearly all sky, and 
if the sky noise behaves as a Gaussian, the slope of 
the cumulative pixel histogram is l/(-\/27rcr) which 
can be used to find a for the individual cells. The 
final a for the plate is then found from the mode 
of the histogram of the individual cell a's. This is 
a good approximation for a over the whole plate 
from the unsmoothed sky using all pixels in the 
plate. Finally, the edges of the average sky image 
are examined in an attempt to avoid processing 



the "bad" edges of the plate as well as the label 
and sensitometer areas: areas in which the mean 
values are more than 5 a away from the mean sky 
are set to a negative value, which is a flag for the 
object detection program to avoid processing. 

3.2.2. Object Detection 

This task used the same implementation of the 
ROE COSMOS image processing software used 
for the GSC-I fseelLasker et al.lll99nl : lLutJll98(i 
MacGillivrav and Stobid 119841 ) which connects a 
minimum number of pixels that are brighter than 
a threshold level above the sky. The threshold 
level for object detection was nominally set at 3a, 
though it was increased to 4ct in cases of plates 
with large grain noise or in crowded fields. A sin- 
gle value of sigma was used for the entire plate, 
as described above. Fifteen additional thresholds, 
equally spaced (in density space) between the de- 
tection threshold and the saturation level, were 
also defined, and minimum areas were determined 
at each of these thresholds as well. A small im- 
age " cutout" was then extracted and saved around 
each detected object, with ample margins to in- 
clude the wings of extended objects. An example 
of image cutouts is given in Figure [H 

From the size, integrated density, and moment 
values of each object, parameters such as object's 
shape and {x, y) of its centroider were derived. 
These and any further image processing were per- 
formed on the cutout images rather than the full 





Fig. 4. — Cutout of a bright star (left panel) and of a resolved binary star (right panel). The irregular 
truncation of the peak of the brightest object is due to the saturation limit of the photographic emulsion. 





Fig. 5. — Typical high latitude field of 4' x 4' (left panel) showing various stellar and extended objects, 
plus a deblended binary system (zoomed in the right panel). The identified objects are encircled by ellipses 
whose size and orientation arc defined by the moments of their density distribution; the line color indicates 
the object classification: green for "stars", yellow for "nonstars", and black for either "defects" or "blends". 



scan image. The maximum cutout extracted was 
4.25' X 4.25', so, if a detected object was larger 
than this, it received no further processing (de- 
blending, centroiding, etc). This results in a small 
number of "holes" in the GSC-II catalog around 
the brightest stars, galaxies, and the occasional 
very bright satellite trails. 

3.2.3. Blend resolution 



Th is is based upon the published (jBeard et al. 
1990l ) ROE COSMOS deblender algorithm, which 
works well also in relatively crowded regions. Each 
new object is re-thresholded, and at each new 
threshold the connec t ivity algorithm is applied 
(JThanisch and RobinI (|l984[ )) to check if the ob- 
ject has split into two or more objects. Faint 
outlying pixels around blended objects are frac- 
tionally allocated according to the faint pixel allo- 
cation technique as described in these references. 
The deblender uses the same thresholds as the ob- 
ject detection task, but applies a single local sky 
value, i.e., the one recorded by the object detec- 
tion. Once an object has been identified as sin- 
gle, its size, integrated density and moments are 
computed and parameters such as x,y and shape 
finally calculated. An example of deblended ob- 
jects formed by the two companions of a parent 
binary is shown in Figure El 

3.2.4. Centroiding 

The centroider algorithm perforins a fit to the 
cutout image with an elliptical Gaussian to find 
the object's center. After the image is normal- 
ized, the centroid is calculated iteratively, with the 
barycentric position as starting point, and assum- 
ing a flat sky. The centroider returns the ampli- 
tude, X, y, and their covariance matrix, as well 
as the X slope, y slope, and level of the fitted 
background. Deblended images occurring in the 
same cutout require special handling. First, an 
attempt is made to perforin a simultaneous fit to 
all objects in the deblended family. If this fails to 
converge, each deblended component is run indi- 
vidually through the single centroider after apply- 
ing the fractional pixel allocation determined by 
the deblender. If the centroider failed to converge 
on an object after a reasonable number of iter- 
ations, or if the centroider failed for some other 
reason, only the barycentric position determined 
from the object detection or dcblending was used 



in the calibration. If the object had been classified 
as a plate artifact, which has a very low chance of 
convergence, the iterative procedure was not at- 
tempted. For objects that occurred in very large 
cutouts (for example, faint stars in the vicinity of 
a bright star), the centroider was also omitted, due 
to resource constraints. 

3.2.5. Classification Parameters 

In order to provide additional information for 
the image classifier (see Table [2|), two additional 
types of numerical features known as spike pa- 
rameters and texture parameters are computed. 
Bright point sources on Schmidt plates display 
diffraction spikes. These were used to good ef- 
fect in GSC-I for star/nonstar discrimination in 
the brighter magnitude ranges. Spike parameters 
result from the comparison of the image gradient 
along the spike, which should be weak, with the 
gradient along the thin dimension of the spike, 
which should be strong. Texture parameters are 
based on the co-occurrence matrix of the object 
cutout (with fractional pixel allocati on applied 
i n the case of deblended objects), see iMalagninil 
(|l983l ) for further details. 



3.2.6. Summary of Known Limitations 

The sky processor is less efficient in very 
crowded fields or around bright stars and halos, 
such as filter ghosts. As already mentioned, the 
object detection procedure fails in the presence 
of large objects, such as bright stars or galaxies, 
or connected objects that will require a cutout 
greater than 4.25'x4.25'. Also, the use of a unique 
sky standard deviation for the whole plate be- 
comes problematic at some low galactic regions 
in the presence of high extinction gradients; here, 
local detection thresholds should be determined 
for optimal object detection. 

GSC-II has both long- and short-exposure 
plates along the galactic plane. In these regions, 
the plate processing pipeline was often run with 
customized parameters. 

Faint objects in the halos of bright ones were of- 
ten missed due to both the sky problem described 
above and a failure of the deblender task. Objects 
brighter than 13th magnitude were sometimes bro- 
ken up into smaller objects. The halo or spikes of 
a bright star were usually split up and identified 



as artifacts; if these artifacts were matched be- 
tween successive plates or if they were classified 
as star/nonstar objects, then they were included 
in the export catalog. In some cases, faint objects 
close to a bright stars have been erroneously de- 
tected twice: firstly, as a secondary component in 
the wide cutout of the bright star, and secondly 
as a direct detection. 

3.3. Plate calibrations 

For each plate used in the GSC-II, the results of 
the image processing pipeline discussed in the pre- 
vious sections provide a list of identified objects 
with position measurements (and errors), inte- 
grated density measurements, and a set of shape- 
related parameters. These are the raw data used 
by the calibration procedures to derive celestial co- 
ordinates, natural magnitudes, and classification. 

3.3.1. Astrometry 

The astrometric calibrations of the GSC 2.3 are 
based on reference catalogs that provide local rep- 
resentations of a fundamental coordinate system. 
The chosen astrometric reference catalogs are the 
ACT and Tycho-2, as described in Section 2.2. 
Since, on average, there are a few hundred ref- 
erence stars per plate, the plate positions are tied 
directly to the reference catalog and hence to the 
ICRS. The brighter stars appearing in the refer- 
ence catalogs, i.e., V < 8.5, are heavily overex- 
posed on the Schmidt plates and thus not useful 
in calibrating the plates. 

In order to compensate for differential refrac- 
tion across the plate, the measured x and y pixel 
coordinates are pre-correc ted for refract ion using 
the method described by iKonigi ( 19621 ). Deter- 
mining the relationship between x, y and their as- 
sociated equatorial coordinates was done by tra- 
ditional global plate modeling, after an equidis- 
tant projection was applied to transform equato- 
rial into tangential coordinates. We found that a 
quadratic model was suitable to represent the rela- 
tionship between the a; and y plate measures (pre- 
corrected for refraction) and the equidistant tan- 
gential coordinates (^, 77). Only the objects flagged 
as astrometric reference stars were used in a least- 
squares reduction to determine the coefficients of 
the model. These coefficients were then used to 
map all the x, y measurements to the associated 



(^, rf) values. Having started calibrating the plates 
in 1998, when Tycho-2 was still not available, the 
polynomial plate model was determined from the 
ACT for w 73% of the plates and from Tycho-2 
for the remaining. 

Once a significant number of plates from the 
same survey had been reduced (approximately 
100), the estimated tangential coordinates were 
compared to those calculated from the reference 
star's coordinates, and used to correct for residual 



system atics according to the precepts of lTaff et al 
( 19901 ). Specifically, for each plate the reference 



stars' residuals were accumulated in 4.4 mm x 4.4 
mm bins, then added and averaged over all the 
plates to build a plate-based astrometric mask, 
which yields the systematic pattern the global 
plate model failed to remove. A separate astro- 
metric mask for each survey was determined using 
the ACT catalog. Figure [6] shows the mask applied 
to the northern blue survey plates (XJ). 




Fig. 6. — XJ astrometric mask 40x40 bins. The 
scale of 1 arcsec is plotted in the lower right-hand 
corner. 



3.3.2. Photometry 

GSC-II magnitudes were derived via modeling 
of the nonlinear density-to-intensity response of 
plate-based photometric calibrators, in the natural 
system defined by each individual plate (Table [J). 
The calibrators are essentially based on Johnson- 
Kron-Cousins B, V, R standards from the Guide 



Star Photom etric Catalog (GSPC) I (jLasker et al. 



19881 ) and II (IBucciarelh et al.ll2001l ): additionally, 
Bt, Vt photometry of Tycho stars was included in 
the calibration to constrain the bright range of the 
response curve. 

Initially, the magnitudes of photometric cal- 
ibrators have been transformed from the origi- 
nal Johnson-Kron-Cousins and Tycho photomet- 
ric systems to the natural system of the plates 
by means of linear interpolation of the color- 
color point diagrams shown in Figure [71 These 
were generated by synthetic photometry (IR AF 
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SYNPHOT package lBushouse and Simonlll994D of 

stellar spectra from the Bruzual-Perss on-Gunn- 
Stryk er Spectrophotometry Atlas (Gunn and Strvker 
1983I ). Note that in the usual case of photometric 
standards having both B-V and V-R available, 
the transformation based on the color closer to 
the requested photographic passband was utilized 
(e.g. Bj — Vpg and T^g — Rp as a function of B — V 
and V — R, respectively). 

The instrumental photometric parameter of 
choice was the integrated density, defined as the 
logarithm of the total density. 



P = loglO Y] {Di.3 - ^bck) 



(1) 



2J 



calculated as the sum of the PDS-measured 
density, Di^j (ADC units), above the local back- 
ground, -Dbck, for each pixel pertaining to the se- 
lected object. The fit to standard magnitudes was 
performed using Chebyshev polynomials, namely: 



4 



TkiP'), 



where To(x) = 0.5 , Ti{x) = x, T2{x) = 2x^ -\, 
T-i{x) = 4a;3 - 3a;, and Tiix) = 8a;'* - 8x^ + 1, 
Ak are the unknowns to be estimated, and P' is 
the photometric parameter linearly scaled to map 
the range [-1,1], where Chebyshev polynomials can 
be suitably applied. A typical pipeline output, 



as given in Figured shows that these polynomial 
bases fit nicely the characteristic curve over a large 
range of densities, in particular at the faint end, 
where the large majority of the objects lie. In this 
respect, extensive experiments have demonstrated 
that the use of polynomials of order higher than 
fourth was not advisable because of their high in- 
stability toward the plate magnitude limit. The 
resulting coefficients were then used to convert the 
measured integrated densities into magnitudes for 
all the objects detected on the plate. 

As a protection against gross errors in the cal- 
ibration, reference stars with large residuals were 
rejected from the fit, according to the criterion A 
> (7 X ei , where A is the residual of the star under 
test, (7 is the formal error of the fit, and ei is an ad- 
justable parameter, generally set equal to 2. The 
procedure was stopped when no further stars were 
rejected, or when the difference of formal error be- 
tween successive iterations was less than £2, an 
adjustable parameter typically set at 0.03. Occa- 
sionally, such iterations failed to achieve formally 
acceptable solutions. In this case, manual inter- 
vention was required, for example, to change the 
order of the Chebyshev polynomial or to modify 
the acceptance criteria for the reference stars. 

The color transformations become poorly de- 
fined at extreme color indices; therefore, only stan- 
dard stars in the ranges of —0.34 < B — V < 1.85 
and -0.14 <V - R< 1.42 were accepted. More- 
over, those stars for which the integrated den- 
sity could not be adequately determined were dis- 
carded from the initial set. This could depend on 
a variety of reasons, such as image processing fail- 
ures, defective images, and centroiding biased by 
faint components of unresolved blends. Therefore, 
in general, stars brighter than 8.5 mag as well as 
deblended objects were eliminated. Since Tycho 
stars are generally more numerous than GSPC-II 
stars on any plate, thereby overconstraining the 
bright side of the fit, a weighting scheme deter- 
mined by the magnitude range of the selected cat- 
alog divided by its number of stars was applied. 
Figure [9] shows an example of the spatial distri- 
bution of the photometric calibrators, with the 
GSPC-II sequence near the plate center and the 
Tycho stars located all over the plate. 

Generally, GSPC (I/II) stars range between 
9th to 19th magnitude (in V), while the GSC 
objects are both brighter and fainter, such that 
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Fig. 7. — Color-color plots used to transform GSPC-I/II calibrators from the BVRc photometric system to 
photographic passbands Bj, Rp, and In (x — y lines are displayed for reference). 
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Fig, 8, — Photometric calibration of plate XP443, Small circles represent GSPC-I/II and Tycho stars used 
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applied in this case. 
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Fig. 9. — Photometric calibration of plate XP443. Plate distribution of the photometric calibrators (left 
panel). Note the faint GSPC-I/II sequence located at X « 8000 and Y w 10,000, while Tycho stars are 
spread over the plate. The right panel shows the objects rejected for various reasons (blends, outliers, colors 
outside the range, etc.), plus a GSPC-I/II sequence very close to the plate border {X w 22, 000, Y sa 18, 000). 
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Fig. 11. — Distribution of total density vs. peak 
density for objects from Figure [10] after transfor- 
mation using ranks. Stars and galaxies from differ- 
ent photographic plates separate very well. An ac- 
curate classifier can be constructed based on these 
ranked parameters. 



extrapolations from the photometric reductions 
are required. Because the images of bright stars 
(brighter than 9.0) are heavily exposed on Schmidt 
plates, their integrated density is not well deter- 
mined. Therefore, no extrapolation of the pho- 
tometric reduction for objects brighter than those 
used in the calibration was performed, and the cor- 
responding magnitude flag set to —99. In the lat- 
est catalog version, all the bright magnitude stars 
having a Tycho-2 counterpart were replaced with 
the Tycho-2 magnitude value (see Section 4). 

The method of extrapolation used for the faint 
stars was tangential, i.e., a linear extrapolation 
with the slopes set to be tangent to the end points 
(brightest and faintest standard stars) of the fitted 
function. However, by using a fourth-order Cheby- 
shev polynomial, the extrapolated region was of- 
ten unstable. In order to stabilize this region, we 
introduced an artificial reference star and assigned 
to it the plate magnitude limit of the survey — 
obtained from the literature and reported in Ta- 
ble[l] — which was associated to the integrated den- 
sity equal to that of the smallest detectable image 
on the plate. In general, the use of an artificial 
point improved the situation, especially on plates 
where the faintest GSPC stars were in the range 
of 16-17 mag. On the other hand, in cases where 
the standard sequence had a fainter limit, the use 
of the artificial star was not needed and could 
even cause GSPC stars to be rejected by the it- 
erative process — an undesirable effect — so the ar- 
tificial point was not used on these plates. 

Variations in the photometric sensitivity across 
the Schmidt plates are caused by sensitivity varia- 
tions in the emulsion, telescope vignetting, and 
point-spread function (PSF) changes across the 
field. We have not corrected for such effects with 
this release. 

It should also be emphasized that the photo- 
metric pipeline, which is tuned for point-like ob- 
jects, systematically overestimates the magnitude 
of bright galaxies {Rp < 18) which do not have a 
stellar PSF. This limitation was acceptable since 
the primary purpose of GSC-II is to provide guide 
stars for telescope operations. It is possible that 
a future version may use a different calibration 
technique to provide more accurate galaxy magni- 
tudes. 
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3.3.3. Image Classification 

The main motivation for classification in the 
GSC series is to reliably identify stars, which can 
be used for guiding, from "nonstars" , which can- 
not. At the scale of the Schmidt plates, such 
"nonstars" include not only galaxies, but blends of 
overlapping images which were not resolved into 
individual objects by the deblending step of the 
plate processing pipeline. A third category, "ar- 
tifact", is used to label images that appear on 
the plates but do not correspond to astronomical 
objects — for example, pieces of diffraction spikes, 
halos, etc. In a subsequent release of the catalog, 
we plan to replace the generic "nonstar" category 
by "galaxy" and "blend" labels, at least for the 
brighter objects where such discrimination is pos- 
sible. In general, nonstars can be considered to be 
primarily galaxies far from the galactic plane, and 
primarily blends close to the plane. 

A set of 30 features, describing the shape, uni- 
formity, and brightness, was calculated for each 
object, as detailed in Table [21 This set was se- 
lected based on good star/galaxy separation in 
feature plots. The "raw" features were then trans- 
formed by the statistical process of " ranking" , that 
is, ordering a population of objects by the value 
of the feature, and mapping the full range of the 
feature into a range of 0:1. The ranked features 
were then used i n the remain der of the classifica- 
tion process fsee IWhitelll997 ). This ranking has 
the effect of significantly reducing plate-to-plate 
variations in the features (see figures fTOl and fTTjl . 
removing the necessity of creating a specific train- 
ing set for each plate. This can be illustrated 
by considering the case of ellipticity, a common 
classification feature: the roundest objects on any 
plate are likely to be stars, even if one plate is 
poorly guided compared to the other. The rank- 
ing was performed separately in different zones of 
the plate, which was divided typically into a 7x7 
grid, whereas some especially crowded plates were 
gridded more finely due to memory limitations. 
Objects lying near the grid boundaries were inter- 
polated. 

An oblique decision tree (OCX iMurthv et al. 



19941 ) was used to perform the classifications. Un- 
like classical, or axis-parallel, decision trees that 
can split on only one variable at a time, OCl can 
perform splits on a linear combination of all the 



features. This permits OCl to find and sensibly 
exploit relationships among the image features. To 
lower the variance contribution to the classifica- 
tion error, a set of five decision trees were created 
from the training set of 5334 hand-classified ob- 
jects based on stars and galaxies from the deep 
CCD catalog of IPostman et al.l (|l996l ). using dif- 
ferent randomizations in the tree-building process 
to produce five pruned trees from the same data. 
During the plate classification task, each object 
was classified by all five trees, and the results were 
voted to produce a final, single plate, classifica- 
tion. 

The training set and decision trees were tuned 
to produce the fewest number of misclassifications 
over the entire population of objects on the plate. 
Since the majority of the detected objects are faint 
objects, and since bright and faint objects clas- 
sify differently, this effectively meant optimizing 
for faint objects at the expense of bright objects. 

4. The GSC-II Database 

The COMPASS database is the primary repos- 
itory for all the data from the plate image pro- 
cessing pipeline. The GSC-II is exported from the 
database as a subset of the multiband and multi- 
epoch surveys because the same objects on the 
sky appear on different plates. The database es- 
tablishes the associations between objects on the 
different surveys and provides a central indexing 
scheme for which the export task can extract the 
calibrated parameters for each object. 

4.1. Database Structure and Operations 

The database architecture is built on a hierar- 
chical object-oriented framework using Objectiv- 
ity, a commercial Object-Oriented Database Man- 
agement System. We regarded as important ele- 
ments for the selection of this system its low devel- 
opment costs and its scalability to multi-terabyte 
data sets, along with the capability of mapping 
the object modeling of the raw data directly into 
the database schema. The latter has allowed us 
to identify a unique correspondence between the 
plate-based raw data and partitions of the celestial 
sphere (sky regions) through the creation of a fed- 
erated database file system (see diagram in Figure 
I12p , with each plate covering multiple regions and 
each region containing several overlapping plate 
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sections. 

The subdivision of the celestial sphere is based 
on the HTM (Hi erarchical Triangulated Mesh 
iKunszt et alj 12001 ). a spherical quad-tree index 
library developed by the Sloan Digital Sky Survey 
(SDSS) project at Johns Hopkins University. The 
regions are approximately equal area of the order 
of 1 square degree. Several large-scale archives 
have adopted HTM to facilitate interoperability 
for data correlations in addition to providing a 
high-speed spatial index into a data set. 

Each of the HTM regions in the database con- 
tains a number of "containers" corresponding to 
the plate measurements, the reference catalog data 
(GSC-I, GSPC-H, Tycho-2, ACT and SKY2000), 
and a master index which identifies each unique 
object on the sky within that region and has links 
to the individual observations and catalog entries. 

4.2. Object Matching and Naming 

As each plate was processed and loaded into the 
database, a matching task was run which identi- 
fied objects on that plate with previously loaded 
objects. Object matching was based on position 
only and was run iteratively with search toler- 
ances increasing from 1 to 4 arcsec. The match- 
ing code computed the distance with respect to all 
the previous observations belonging to any source 
within the searching radius and assigned the new 
object to the source containing the nearest neigh- 
bor. This procedure is repeated with a larger radii 
on the remaining objects and unmatched objects 
are added as new sources to the master index. 

A known problem with the matching is due to 
noise fluctuations in the extended saturated cores 
of very bright stars resulting in the detection of 
spurious objects. Typically, these bright stars are 
contained in the SKY2000 and Tycho-2 catalogs, 
so the matching of these spurious objects to the 
same catalog entries is the cause of multiple entries 
in the GSC-II. This is not a problem for HST op- 
erations since the r2/c/io-^/SKY2GGG values take 
precedence but we plan on fixing this issue in a 
later release. 

In order to maintain backward compatibility for 
HST operations, the master index was preloaded 
with the GSC-I objects so that the GSC-I identi- 
fications were maintained. As objects were added 
to the index, an internal identification was as- 



signed which was the HTM container name (e.g., 
N012301) plus a running sequential number. This 
was adopted to aid fast object lookups within the 
COMPASS database. It should be noted that this 
name is converted to a 10-character base-36 num- 
ber known as the HSTid in GSC2.3 and later and 
that this is the "ofhcial" name. This was done be- 
cause of the ten-character limit for the Guide Star 
names throughout the HST ground systems and 
flight software once HST pointing switched over 
to using the new catalog for pointing. 

4.3. The Export Catalog 

The GSC 2.3 export catalog provides positions, 
multicolor photometry, and morphological classi- 
fication for about one billion unique objects. For 
each source, the GSC-II database may contain 
multiple measurements deriving from the obser- 
vations matched from the overlapping multi-epoch 
and multicolor plates. In this section, the selection 
criteria adopted to filter the available data are il- 
lustrated. It is important to recall in this context 
that the main driver for the GSC-II project was 
to support and/or to improve on the efficiency of 
observation planning and operations in space and 
ground observatories which sponsored the project; 
the catalog extraction rules ultimately define the 
properties of the catalog released. 

The magnitude limits of the export catalog were 
set at Bj <19.5 or Rp <18.5 for the first public 
version, GSC 2.2, released on 2001 July, while for 
the GSC 2.3 almost all the objects down to the 
magnitude limits of the plates (see Table [1]) have 
been extracted from the database and included in 
the export catalog. 

The magnitude limits of the GSC 2.2 basically 
correspond to the range which satisfies our guide 
star requirements for telescope operations in terms 
of positional and photometric accuracy. In partic- 
ular, the colors provided in the export catalog are 
indispensable for both bright-object protection of 
the MAMA detectors on HST and adaptive op- 
tics operations of ground-based telescopes. The 
classification between stellar sources and nonstel- 
lar objects (extended objects or blends), is critical 
for the selection of both guide stars and adaptive 
optics reference stars and the GSC-II is sufficiently 
reliable within the stated magnitude range. 

The deeper GSC 2.3 release was mainly mo- 
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tivated by our commitnient to distribute to the 
international astronomical community a large and 
valuable data set for a wide variety of astrophys- 
ical studies. Moreover, GSC-II data have been 
already used for various technical analysis which 
required a complete deep all-sky catalog, such as 
the guide star studie s for the JWST (|Spagnall200ll : 
Stvs and Krisal2003r). and the d efinition studies of 
the Gaia mission (JDrimmel et al .. 2005.). 



4-. 3.1. Source selection. 

All the database entries were screened, but only 
objects detected and measured at least on one II- 
laF or IllaJ or IV-N plate were exported. Ob- 
jects classed as defect on only one plate were dis- 
carded. A magnitude threshold, m < 25, signifi- 
cantly fainter than the plate limits of all the pho- 
tographic surveys was applied in order to prune 
both faint false positive detections and spurious 
objects. 

In order to improve the completeness and accu- 
racy of the bright guide stars, which are strongly 
saturated in the GSC-II long-exposure plates, the 
GSC 2.3 has been supplemen ted with data from 
the Tycho-2 and SKY2000 (M vers et all2002f) cat- 
alogs which contain stars down to V ^ 12 and 8, 
respectively. The data for any given object were 
taken preferentially from the Tycho-2 when avail- 
able, then SKY2000 and finally from the Schmidt 
plate measurements. 

The GSC 2.3 also contains the names of the 
GSC-I objects that have been matched with the 
GSC-II objects. 

4.3.2. Astrometry. 

Generally, GSC 2.3 provides positions, {a,S), 
measured on the red IIIaF plate, and multiple 
measurements are resolved by selecting the closest 
solution to the plate center. In case the object was 
not detected on any IIIaF plate, then the positions 
from the blue Ilia J, visual IlaD, and infrared (IR) 
IV-N plates have been used following this priority 
order. The position epoch is set accordingly. The 
position errors are estimated by means of a sim- 
ple general model aimed to provide conservative 
uncertainties for telescope operations. These fig- 
ures do not correspond to formal errors and cannot 
be used to define confidence limits in a statistical 
sense. 



If a source was selected from the Tycho-2 cat- 
alog , the GSC 2.3 provides the positions, the 
proper motions and the mean observing epoch in 
right ascension (R.A.), which is slightly different 
from the declination (decl.) epoch in the Tycho- 
2 catalog. For SKY2000 objects, coordinates and 
proper motions are provided at epoch J2000.0. It 
should be noted that at this time we are not pub- 
lishing computed motions due to a small system- 
atic error in the coordinates (see Section 5.3) 

4.3.3. Photometry. 

Photographic magnitudes Bj, Rp, Vpg, /at, plus 
the blue POSS-I O are provided, if available, for 
all the exported objects. In case of multiple mea- 
surements, which may occur at the plate borders 
where different plates of the same band overlap, 
we selected the observation closest to the plate 
center, without applying any averaging procedure. 
As a result of the main source selection criterion, 
at least one Bj, Rp, or I^ magnitude is always 
present, while additional blue O and Fpg magni- 
tudes are given for the objects also matched on the 
POSS-I O, Palomar Quick-V, and SERC Quick-V 
surveys. As for the positions, magnitude errors 
are approximate and conservative estimates. 

In the case of Tycho-2 objects Bp and Vp mag- 
nitudes are provided, whilst SKY2000 objects may 
provide additional Johnson UBV magnitudes. 



4.3.4. Source Classification. 

Classification voting was implemented on the 
set of observations for an individual sky source. 
For each source we have at least two classifications 
coming from the plates in the different bands, but 
we may have up to 25 if the object falls into a 
region of multiple plate overlap. The final pub- 
lished classification of a matched object is the 
mode of the individual classifications (ignoring ar- 
tifact classification on the assumption that any ob- 
ject matched on two plates is real). Unmatched ar- 
tifacts are excluded from the catalog. In the event 
of a tie, e.g. as many stellar as nonstellar votes, 
we decided in favor of non-stellar based on our ex- 
perience that it is easier to misclassify a galaxy as 
a star rather than vice versa. Individual classifi- 
cations based on older 25 /xm data are excluded 
from the voting to ensure consistency. 
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Each source was classed as star or nonstar and, 
in the case of ties, nonstar classifications were se- 
lected to ensure highest quality stellar sources. We 
used only the classification of the objects derived 
from the higher-resolution scans (1.0" pixel); clas- 
sification of low resolutions scans (1.7" pixel) were 
utilized in the rare case of objects without any de- 
tection on the high-resolution scans. A quality 
flag which indicates an unanimity voting was also 
added to the export catalog. "Defects" revealed 
and matched on more than one plate were classi- 
fied as nonstars. 

Since this catalog is primarily for the selection 
of guide stars and the GSC-I classification was 
strongly biased toward this, then if an object was 
cross-matched with the GSC-I that classification 
was adopted and it superseded the voting scheme 
described here. 

Finally, the morphological parameters (elliptic- 
ity and orientation angle) of the image used to 
derive the positions were exported too. 

4.3.5. GSC 2. 3 Export file format 

The export catalog is a distributed FITS file 
system based on the linear quad-tree HTM sky 
partitioning scheme, with a level-6 trian gulation 
implementation (see iKunszt et al.l 120001 ). Each 
FITS file includes all calibrated celestial sources 
for a single triangular section of sky (HTM level 
6 leafnode), roughly corresponding to one square 
degree of sky data. The file names are the quad- 
tree representation of each triangular section, with 
the first character being either an "N" or "S" rep- 
resenting north or south. 

Each FITS file contains three sections, the re- 
quired primary header, an ASCII table extension 
(containing a lookup table that access software can 
use to speed up reading data subsets), and a BI- 
NARY table containing the actual catalog data. 

Table |3] shows the fields officially released with 
GSC 2.3, including object names, astrometric and 
photometric parameters, object classification, as 
well as some auxiliary parameters and flags. The 
binary FITS tables also include several empty 
fields for future releases or containing data for in- 
ternal testing; these are not part of GSC 2.3 and 
must not be used nor distributed. 

As explained in the notes of Table [H all the 
photometric fields may contain magnitudes for dif- 



ferent, although similar, passbands. GSC-II users 
can identify the exact filter by looking at the 
magnitude-code fields, for which we list in Table 
Uthc photometric codes used in GSC 2.3. 

Further details on the processing of the GSC-II 
objects are given in the source status flag, whose 
structure is described in detail in Table El Tycho- 
2 stars are identified by a Quality Status Flag of 
999999nrii| and SKY2000 stars are identified by a 
Quahty Status Flag of 88888800. 

4.4. GSC 2.3 distribution 

The export catalog, GSC 2.3, requires a stor- 
age capability of approximately 170 GB and it 
is currently maintained by STScI and INAF- 
OATo. Copies of the catalog have been already 
distributed to the project patrons (ESQ, GEM- 
INI, ESA) as well as the main astronomical data 
centers and several research institutions. Due to 
limited resources, the GSC-II team does not plan 
to provide an extensive distribution service of the 
complete catalog and requests will be reviewed 
on a case- by-case basis. However, public access 
to GSC-II is available through both Web and 
VO-compliant (Virtual Observatory) interfaces at 
STScI, OATo and CDS. 

5. Catalog Properties 

The application of the "selection" criteria ex- 
plained in the previous section generated the ex- 
port catalog, GSC 2.3, containing 945,592,683 ob- 
jects in total, classified as stars (22%) and nonstars 
(78%). 

Tableinipresents the global statistics of GSC 2.3, 
including the number of objects for which photo- 
graphic magnitudes, Bj, V^g, Rp, In, are pro- 
vided, and the number of objects related to GSC-I, 
Tycho-2, and SKY2000. 

All-sky cumulative starcounts in the Bj, Rp, 
In are provided in Table [T) 

The GSC2.3 all-sky density map in galactic co- 
ordinates is presented in Figure fT3l which includes 
both stellar and nonstellar objects. The counts are 
dominated by stars close to the galactic plane, and 
most of the objects classified as nonstars are actu- 
ally unresolved stellar blends belonging to crowded 



^Here nn is the number of observations used to compute 
Tycho 2 position and proper motion. 
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fields toward low galactic latitudes. Also, an anal- 
ysis of the counts shows that at Rf <20 about 
4% of the GSC2.3 sky has a density of above 10^ 
objects per square degree, reaching in some cases 
up to 10^ before getting confusion Hmited. 

The reliability of the astrometry, photometry, 
and classification of GSC 2.3 can be investigated 
by comparisons with independent reference mate- 
rial as well as by internal checks. The latter take 
advantage of the significant plate overlap and can 
probe the degree of catalog consistency down to 
the very plate limit; nonetheless, they are not dis- 
cussed here because, since the plate-based nature 
of GSC-II is not preserved in the way the export 
catalog was constructed, an error analysis based 
on the plate overlaps would not give the right in- 
formation about the GSC 2.3 errors. 

5.1. Comparison Catalogs 



The UCAC 2 (jZacharias et al.l 120041 ) contains 
almost 50 million objects down to i? « 16, cov- 
ering the entire southern hemisphere and a large 
fraction of the northern one (complete to 5 = 40°, 
and up to 5 = 52° in some areas). Its astrometric 
precision is between 40 and 70 mas per coordi- 
nate depending on magnitude, and systematic er- 
rors are within 10 mas; these values are lower than 
the typical GSC-II errors, except for the faintest 
UCAC 2 objects, so that this catalog appears ideal 
to reveal the astrometric precision of GSC 2.3. 
The only drawback is the magnitude limit being 
significantly brighter than that of GSC 2.3. It is 
worth noting that the next release of this catalog 
(UCAC 3) which is all-sky would be an excellent 
reference catalog to recalibrate the plates since it 
is deeper than Tycho-2 and would reduce some of 
the systematic errors in GSC-II that we describe 
later. 

In order to extend the astrometric error analy- 
sis to fainter mag nitudes we also matched GSC 2. 3 
with SDSS DR5 (|Adelman-McCarthv et all2007l) . 
This provides astrometrically and photometrically 
calibrated data for an area of ^ 8000 deg^ around 
the northern Galactic cap; a total of about 6 • 10^ 
objects including stars, galaxies and quasars are in 
common with GSC 2.3. For point sources brighter 
than r — 20, the astrometric accuracy of SDSS 
is 75 mas per coordinate with an addition al 20-30 
mas systematic error (see lPier et al.ll2003l ). At the 
survey limit (r = 22), the astrometric accuracy is 



limited by photon statistics to approximately 100 
mas root mean square (rms) . The photometric cal- 
ibration is accur ate to roug hly 0.02 mag in the g-, 
r-, and i-bands (jlvezic et al.ll2004l) . so the SDSS 
can also be used to estimate the GSC 2.3 photo- 
metric errors. A comparison against GSPC-II also 
provides a final global validation of the photome- 
try. 

In order to compare the astrometry to an 
all-sky catalog of similar size, yet derived from 
independent data, we also matched GSC 2.3 
against the Two Micro n All Sky Survey (2MASS; 
ISkrutskie et"aD (|2006h ) Point Source Catalog 
(PSC). This lists 470,992,970 objects distributed 
over the entire sky, providing near-IR photome- 
try with signal-to-noise ratio (S/N) better that 
10 down to J ~ 14.2, H ~ 15.1, and Ks ~ 14.2. 
Note that the GSC 2.3 does include the bulk of 
2MASS objects, as these basically span the entire 
magnitude range of GSC-II, thanks to the wide 
color spectrum covered by both catalogs and the 
natural spread of stellar temperatures and galactic 
extinction. The astrometric accuracy of 2MASS 
is 70-80 mas per coordinate over the magnitude 
range 9 < A's < 14, while for brighter sources 
it is approximately 120 mas. For fainter sources 
the error increases monotonically. The systematic 
errors are of the order of 10 mas, on average, and 
up to a maximum of 25 mas for the worse cases. 

The USNO-B (JMonet et al. l l2003h or the South- 
ern SuperCosmos Survey (SSS, iHamblv et al.l 



20011 ) were not used as external checks since they 



are derived from much of the same plate mate- 
rial as GSC-II, so do not provide independent 
measurements of the quality. As a sanity check 
GSC 2.3 was compared to both USNO-B and SSS 
but the differences did not reveal anything new 
that is considered significant. 

In addition to the large-scale comparisons 
performed with the above catalogs, we further 
checked object classification against the selected- 
by-eye galaxy c atalog of t he NPM/SPM programs 
(|Klemola et al-l flggi: Pla tais et al.ll2002l) and var 



ious NGP surveys (llnfante 
the APM (Irwin 



2007) 



^993; Gr aham et al. 
and McMahonI 1 19921 : 



Lovedavl [l996l ) and APS (|Cornuelle et al.l [1991 



automated galaxy classification catalogs. 
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Fig. 13. — GSC 2.3 all-sky map. Cumulative counts in galactic coordinates, including both stellar and 
extended objects. The color scale indicates the GSC 2.3 density ranging from to 60,000 objects per square 
degree. The image resolution is a smooth version of the densities obtained from the HTM counts. 
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Fig. 14. — Comparison to the UCAC 2 catalog for a sky region centered at h, spanning all declinations, 
and in the magnitude range 14.5-15.0. The panels show the results of the x^ (top panels, the histogram is 
the theoretical x^ distribution), and K-S test (bottom panels, the dotted line is the theoretical distribution). 
The empirical distributions do not satisfy the tests (S/Q > 1, see the text). This is attributed to the presence 
of residual systematic errors. 



5.2. Astro metric errors 

As discussed in Sect. 3.3 (astrometric calibra- 
tions), the GSC-II is affected by systematic errors, 
which are function of both position and magni- 
tude. In the following, we investigate size and vari- 
ation of such errors by comparing the GSC2.3 as- 
trometry to the positions of selected external cat- 
alogs. These were calibrated with Tycho-2 stars, 
thus we are in a position to check at the same 
time on the degree of global and local adherence 
of the GSC 2.3 positional system to that of the 
ICRF reference frame. We also give an estimation 
of the positional random error for both point-like 
and extended objects, which is the relevant quan- 
tity in any application requiring knowledge of the 
relative astrometric precision. 



5.2.1. Statistical tests 

A purely statistical approach to test the de- 
gree of randomness of the astrometric errors is the 
method of hypothesis testing. Specifically, we ap- 
plied on our data two goodness-of-fit tests, i.e., 
the x^ and the Kolmogorov-Smirnov (K-S) tests, 
as detailed in the following. 

Let €ii and 6^2 be the actual coordinate er- 
rors of star i from the GSC 2.3 and the compari- 
son catalogs, respectively. Assuming that en and 
ei2 are independent normally distributed random 
variables with respective variance ctJ and cr|, then 
their difference di — en ~ ei2 is also normally dis- 
tributed, with variance equal to a^ + crl- 

This assumption on the distribution function of 
di was the working hypothesis to which the x^ and 
the Kolmogorov-Smirnov tests were applied. Since 
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Fig. 15. — As in FigurJTil these panels show the results of the x^ and K-S tests for the GSC 2.3 versus 
UCAC 2 residuals, this time restricted to one plate and to the magnitude range Rf — 14.5 — 15. In this 
case, both significance tests are satisfied [S/Q < 1). 



each star has only one observation per plate, we 
had to resort to an "idealized" experiment in order 
to have the required redundancy. To this end, we 
considered measurements of different stars close in 
magnitude as belonging to the same parent popu- 
lation, i.e., each measure was treated as a single re- 
alization of the statistic under test. Such approach 
is correct as long as the measurement/model er- 
rors are independent (or marginally dependent) on 
sky/plate position or other effects not taken into 
account by the plate model; if any of these condi- 
tions is not satisfied, the goodness-of-fit tests are 
doomed to fail, whereas a successful test would 
confirm the existence of such conditions. 

Assuming a Gaussian probability density dis- 
tribution (pdf) of the di, we estimated first the 
population standard deviation by using the for- 
mula a ~ 



^d, where d 



^J2^\d,-{d)\ is the 



N 



sample mean deviation, which is less sensitive to 



outliers. Then, we applied a 3cr rejection criteria 
to the data and calculated the usual sample mean 
fi and standard deviation a of the cleaned sample. 
Finally, by binning the residuals at 0."05 steps, we 
built the test statistic 



E 



Npo^ 



where N is the total number of residuals, Ni is 
the number of residuals in the i-th bin, and 



(.a:i+Ax/2 

PQi = l/(V27rcr) / exp( 

Jxi~Ax/2 



2a2 



)dt 



the probability of a residual falling into the ith 
of the k bins, which can be evaluated using the 
error function. Note that the a and fi parameters 
used in the integral above are estimated from the 
sample itself. 
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The approximate distribution of x^, valid for 
large samples, is that of the chi-square probability 
function with fc — 3 degrees of freedom. 

An alternative significance test we have applied 
to our sample is the double-sided K-S statistic, 
which is exact in the sense that does not require 
data binning. 

The K-S test statistic is defined as 



Dk = max\Sk{x) — F{x)\ 



with 



Skix) 



{0 for X < di 

k/N for dk < X < dk+i 
1 for dk < X 



and 



F{x) = l/(\/27rcr) 



exp( r-^ )ds 



2f72 



0.5 



erf(a;) 



where erf is the error function. 

Both tests adopted a 5% risk level of rejecting 
the null hypothesis {Hq) of Gaussian distribution 
of the residuals, and therefore compared the statis- 
tics derived from the sample (S) with 95% quan- 
tiles {Q) taken from the x^ with fc — 3 degrees of 
freedom and the K-S theoretical functions. If the 
sample values of the two statistics, 5^2 or Sd„ , 
were larger than the quantiles, i.e. S/Q > 1, the 
Ho hypothesis was rejected; vice versa, if S/Q < 1, 
the Gaussian test was accepted. 

We have tested the residual differences com- 
puted against SDSS and UCAC 2 inside R.A. 
strips of 2 h and spanning the entire decl. range. 
The results of the tests for different magnitude 
bins and for a number of different plates show 
that, although not strictly Gaussian, the distri- 
butions are well behaved. As a general trend, the 
errors in decl. are better behaved, and the R.A. 
error increases more toward fainter magnitudes. 
As an example, Figure [14] shows the results ob- 
tained by comparing GSC 2.3 coordinates against 
UCAC 2 on a strip of sky between and 2 h and 
-90 < 5 < +50 degrees. 

If we do not bin in magnitude, the distribution 
does not behave like a Gaussian as would be ex- 
pected from combining different Gaussian popula- 
tions; however, even limiting the magnitude range 



to within Rp — 14.5 — 15, which gives a total num- 
ber of matched (stellar) objects of about 44,000, 
both tests still fail. In particular, the distribution 
of AacosJ diverges from the expected one, and 
presents a significantly larger dispersion than the 
A5 distribution. This occurrence is indirect ev- 
idence of local systematics errors which distorts 
the shape of the estimated distribution function. 
In fact, if we restrict the sample to a single plate, 
and the same magnitude range (862 stars sample), 
both tests are better behaved, and the Hq hypoth- 
esis is accepted, as illustrated in Figure fTSl 

5.2.2. Magnitude-Dependent errors 

Morrison etHI (|l996l . l200ll ) found that posi- 



tional differences of GSC-I objects matched on 
overlapping plate areas lead to the presence of 
a radial magnitude term with the origin in the 
plate center. Such an effect was characterized 
and removed fror n GSC-I by means of the Astro- 
graphic Catalog (jUrban et al.lll998[) . whose mag- 
nitude limit is around 14-14.5, hence suitable to 
probe the faint limit of GSC-I. However, applying 
the same precepts to GSC-II data would fail to 
detect nonradial magnitude terms which we know 
to exist. 

Monet et al.l ( 2003[ ) were able to use a special 
catalog (YS4.0, magnitude limit V ~ 18) compiled 
from the NPM and SPM plate scans to detect and 
remove magnitude-dependent residual errors dur- 
ing the constr uction of USNO- B by means of the 
mask method (JTaff et al.lll990t ). 

For our magnitude error analysis, we cross- 
matched GSC 2.3 against the adopted external 
catalogs using a search radius of 3 arcsec to recover 
the common objects and build the differences: 



Aa = (acSC - "cat) cos^gsc A(5 = (^gsc - <5cat) 

Then, we binned the data at 0.5-mag steps and 
for each of such bins computed the statistics de- 
scribed below. 

We began with estimating the global rms of the 
residuals as 



4. = ^E(A«i) 



(2) 



where N is the number of objects (stellar or 
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Fig. 16. — GSC 2.3 vs. UCAC 2. Astrometric residuals of GSC 2.3 sources classified as point-like objects 
{left panels) and extended objects (right panels). The symmetric red solid lines show the random errors, 
±(TAa, icTAS, as a function of the magnitude (see Equation [5]) . Asterisks connected by the black dotted line 
show the means (Aa), (AS) averaged over all the HTM regions (Equation [7| , while the error bars are the 
standard deviation of the systematic field-to-field variations (Equation [6|). The orange and blue dotted lines 
show the averaged (Aa), {A6) means over the northern and southern regions, respectively. 



nonstellar) inside a particular magnitude bin. It 
is evident from this definition that eaq represents 
a catalog average of the total error, as it includes 
both systematic and random errors, and it is the 
quantity closest to the real astrometric error. Be- 
sides, relative astrometry over small fields (of the 
order of half a degree or less) is mainly affected by 
random errors; therefore, next we estimated the 
magnitude dependence of the random component 
of the GSC 2.3 error. For each of the 32768 HTM 
regions we measured the local systematics, (Aa)r, 
and its variance, a'^^ ^ , as 



(Aa), = 



-YAc 



Nr 



i=l 



and 



2 



^ (Aa, - (Aa).)' 



(3) 



(4) 



respectively, with Nr being the total number 



of objects per 0.5 mag bin in region r. Conse- 
quently, we estimated the catalog random error as 
the weighted mean of the standard deviations in 
each HTM region, i.e.. 



„2 






(5) 



r = l 



where A^a. is the total number of HTM re- 
gions. Outliers due to mismatches were rejected, 
and GSC 2.3 entries flagged as Tycho-2 stars ex- 
cluded. Also, magnitude bins with poor samples 
(less than 20 matched objects) or very crowded 
samples (more that 5000 objects) were not used 
in this analysis. Analogous quantities were evalu- 
ated for the residuals in decl. 

The statistics in Equations [2] and [5] are sum- 
marized in Table [3 which reports the results of 
the comparisons with UCAC 2 for the bright- 
to-intermediate magnitude range and with SDSS 
for the intermediate-to-faint range. We note that 



24 



stellar Obiects 



Systematic <Aacos6> 
Random a 



fe ttfJH-M +ttff 



Global 
<Aa> North 



14 16 IB 20 

R (mao) 
stellar Objects 



Systematic <A6> 
Random a 



Global 

<A6> North 

South 

14 16 IB ao 

R (ina;^) 





Extended Objects 











X Systematic <Aacos(5> 
Random a^^ 


' 


0,5 




^ : 




f^jl^^f^WW-f+ii : 






0.5 
1 


Global 

<Aa> North 

South - 



R, (mag) 
Extended Objects 




Fig. 17. — Similar to Figure [16] but for the comparison GSC 2.3 versus SDSS. Astrometric residuals as 
a function of the magnitude of GSC 2.3 sources classified as point-like objects (left panels) and extended 
objects (right panels). Same symbols as in Figure [TBI Random position errors (red solid lines) are almost 
constant; for stellar objects ctaq — o^AS — 0.14 down to Rp — 18.5 and comparable with the small scale 
field-to-field systematic zero-point variations shown by the error bars. The dotted lines confirm the presence 
of the systematic offset between the northern/sourthern hemispheres already found at brighter magnitudes 
by the comparison against UCAC 2. Notice that the global mean residuals (asterisks) are very close to 
northern zero-points because of the small fraction of SDSS objects with negative declinations. 
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the SDSS astrometry appears more precise than 
2MASS at the faint limit, and therefore more rep- 
resentative of the actual GSC 2.3 error, although 
for only a portion of the sky. 

The data confirm the better behavior of the 
point-like objects ("stars"), which attain a ran- 
dom error of ctaq — ctAiJ < 0.14 arcsec at inter- 
mediate magnitude, 14 < Rp < 18.5, increasing 
up to 0.23 arcsec at the faint end {Rp ~ 20). As 
expected, larger residual errors are obtained for 
extended objects ("nonstar"), which may include 
galaxies and nebulae, as well as many blends and 
unresolved binaries, in particular at bright mag- 
nitudes or in crowded fields toward low galactic 
latitudes. In this respect, we note that the higher 
errors of the extended objects in the UCAC 2 com- 
parison, as opposed to the SDSS comparison, are 
explained by the dominance of blended images in 
the former sample, whereas the latter is mostly 
made of truly extended objects, being the SDSS 
catalog confined to high galactic latitudes. 

The total rms errors, eAa and e/^s, which also 
include the contribution of systematic errors □, are 
of about 0.2" per component at intermediate mag- 
nitudes (14 < Rp < 18.5), and increase to 0.3" 
down to Rp « 20. 

An estimate of the systematic field-to-field vari- 
ations in each magnitude bin can be obtained by 
taking the standard deviation of the (Aa)r around 
their averages computed over all the HTM regions, 
i.e. 



>{Aa) 



N 

E 



r=l 



A ((A«),, - {Aa)y 



[Nr - 1) 



(6) 



where (Aa) is given by 



r—1 



(7) 



Similar expressions hold for the field-to-field 
variations in decl. The quantities ct/Aq) and cf/as) 
are showed as the error bars in Figures \T^ and [T7] 
and they are close in size to the random errors 
for stellar objects in the intermediate magnitude 
range. The corresponding averages, i.e. Equation 
[7] for R.A., are indicated by asterisks in the same 



^It is convenient to recall here that we chose not to remove 
the errors of the reference catalogs described in Section 15.11 



Figures and represent a uniform weight mean. 
This is a better indicator of the GSC2.3 astro- 
metric zero-point when compared to the straight 
mean taken over all of the residuals of the matched 
objects. For, the simple average tends to overrep- 
resent the regions close to the galactic plane. The 
catalog zero point appears quite small for bright 
objects, as they are close to the magnitudes of the 
Tycho-2 stars used to calibrate the plates, while it 
increases up to ^ 150 mas in R.A. at 20th magni- 
tude. However, a realistic estimate of the GSC 2.3 
adherence to the ICRF must also account for lo- 
cal systematics, which can vary significantly from 
region to region; one large-scale example is the 
northern/southern dicotomy revealed by the or- 
ange and blue lines of Figures [16] and [171 which 
depict the magnitude dependence of the mean as- 
trometric zero point in the northern and southern 
regions. This effect will be addressed in the next 
section. 

5.2.3. Spatial Dependent Errors 

To investigate spatial-dependent errors we have 
checked GSC 2.3 residuals versus SDSS positions 
of 17,507 stars along an equatorial strip with 
145° < a < 205° and for the magnitude range of 
16.5 < Rp < 17.0. From the plots of FigurefTSlone 
can immediately note the plate-to-plate disconti- 
nuity and the gradients within each plate. In par- 
ticular, the systematic zero-point scatter amounts 
to 0.07" and 0.06" in R.A. and decl. respectively, 
as derived by the standard deviation of the run- 
ning mean of the residuals {red solid line), while 
the standard deviation of the residuals is 0.15" and 
0.14". 

An all-sky view of the GSC 2.3 systematic and 
random errors estimated by Equations [3] and |4l 
and utilizing 2MASS as comparison catalog, is 
given by the equatorial maps in Figure [19] for 
the magnitude range of 15 < Rp < 15.5. The 
random error map shows a nonuniform distribu- 
tion, with significant variations with respect to 
the mean values (uAa and cta^) reported in Ta- 
ble [8] In particular, we note the poorer behavior 
along the galactic equator (crowding), as well as 
the dependence on R.A., which might indicate ef- 
fects due to poorer observing conditions (seasonal 
effects and/or suboptimal hour angles during ex- 
posures). The higher errors at the plate corners 
form the peculiar speckled pattern which is most 
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Fig. 18.— GSC 2.3 vs. SDSS. Astrometric residuals (dots), AacosS and A6, for 17,507 stars with 16.5 < 
Rf < 17.0 along an equatorial strip in the range 145° < a < 205°. The red solid lines are the running means 
of the astrometric residuals which indicate the geometric systematic errors. Here, the crosses indicate the 
R.A. of the centers of the POSS-II fields along the celestial equator at i5 = 0. For this sample, the rms's of 
the residuals are 0.15" and 0.14" for R.A. and decl. respectively, while the zero-point scatter is 0.07" and 
0.06" as derived by the standard deviation of the running mean. 
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GSC2.3 vs. 2MASS. a (15.0<R <15.5) - Equatorial coord. 
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GSC2.3 vs. 2MASS. Aa (15.0<R <15.5) - Equatorial coord. 




Fig. 19. — Equatorial maps of the GSC 2.3 regional random errors aAa,r (top panel) and systematic errors 
(Aa) (bottom panel) versus 2MASS, in the magnitude range 15 < Rp < 15.5. The black solid lines trace 
the Galactic equator, while the asterisks represent the Galactic poles and the Galactic center. Similar maps 
are found for the decl. errors. 
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Fig. 20.— GSC 2.3 - 2MASS R.A. (left) and decl. (right) differences as a function of decl. for the Rp 
magnitude ranges 11-13 (top panels) and 18-19 (bottom panels) . The central solid line marks the median 
for the data binned in decl., while the external ones trace the rnis deviations about each median. 
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Fig. 22. — Same measurements as in Figure [^ 
but plotted as a function of decl. to show the sys- 
tematic proper motion error in the southern hemi- 
sphere 



visible along the equator where the effects of the 
sphere mapping geometry are minimized. 

The bottom panel shows a definite systematic 
difference between north and south, which we fur- 
ther analyzed. To this end, we generated plate- 
based vector plots of the R.A. and declination dif- 
ferences versus 2MASS with a 1-mag step reso- 
lution and looked for possible signatures of dif- 
ferent surveys, emulsion/filter combinations, and 
decl. These maps present a characteristic pattern, 
whose typical amplitude is negligible in the mag- 
nitude range of the reference catalog used for the 
astrometric solution, and increases toward fainter 
stars up to an average of ~ 0.25 arcsec (at mag ~ 
20). Moreover, the net residual is positive or nega- 
tive according to decl. /survey telescope, as Figure 
1201 clearly shows. This sign dichotomy, which is 
more pronounced in R.A. than in deck, can be ex- 
plained by noticing that the northeast are flipped 
for a northern versus southern telescope mount- 
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Fig. 21. — Measured proper motions of SDSS quasars in (a) R.A. and (b) decl. Since QSOs should have zero 
proper motion, this is a measure of the proper motion errors (see text). 



ing, while, given the similarities of the observing 
systems utilized (telescopes, filters, plates, strate- 
gies), the systematic components have similar be- 
havior in plate-based coordinates. Therefore, the 
net average astrometric residual as a function of 
decl. changes in sign around the equator for both 
coordinates. As expected, this systematic feature 
becomes more pronounced for fainter stars, far- 
thest from the magnitude range of the reference 
catalog to which GSC 2.3 has been tied. 

The analysis presented in above sections has 
shown that we are dealing with a complex sys- 
tematic pattern of astrometric residuals, contain- 
ing magnitude and spatial dependences ultimately 
leading down to the local nature of the observa- 
tion/reduction process. It is also clear that such 
systematics need to be modeled locally in order to 
be effectively removed; this task will be addressed 
in the next GSC-II release. Notwithstanding, the 
size of the GSC 2.3 astrometric errors is such that, 
although they cannot be neglected for proper mo- 
tion determination, they do not impair the opera- 
tional use of the catalog. 

5.3. Proper Motion errors 

As mentioned in the previous section, we ex- 
pected that the complex astrometric residuals 
would affect production of absolute proper mo- 
tions when derived from simple fitting to the 
multi-epoch positions available in the GSC-II 
database. In order to test this we took a sample of 
quasi-stellar objects (QSOs) that had been iden- 
tified in the SDSS and computed a proper motion 
from the GSC 2.3 measurements at all available 



epochs using a simple least-squares linear fit. The 
QSOs should have zero proper motion since the 
ICRS reference frame is tied to the extragalactic 
dynamical system. The measured proper motions 
are therefore an estimate of our errors. These 
results are shown in Figure [HJ At first glance 
it would appear that the proper motions are de- 
termined with an rms error of 6-12 mas yr^^ (as 
a function of magnitude). However, if one looks 
at these as a function of decl. in Figure [221 we 
see that, whilst the northern hemisphere proper 
motions are determined to 6-8 mas/yr, there is a 
significant systematic error for the southern hemi- 
sphere which is increasing the overall error. We 
believe that whilst the systematic positional errors 
tend to cancel out among the different northern 
surveys, they do not for the southern ones. There- 
fore, although we have exported GSC2.3 proper 
motions from the database, these are used inter- 
nally for telescope operations, and they will not be 
released to the scientific community until the sys- 
tematic component is dealt with and minimized. 

5.4. Photometric errors 

GSC 2.3 photographic magnitudes result from 
the photometric plate calibrations described in 
Sect. 3.3.2, and from the selection criteria which 
defined the export catalog, as discussed in Sect. 
4. Note that the GSC-II pipeline has been tuned 
for point-like sources and it is not optimized for 
galaxy photometry. Similarly to the astrometry, 
GSC-II photometry is affected by both position- 
and magnitude-dependent systematic errors. In 
this section, we assess the GSC 2.3 photometric ac- 
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Fig. 23. — Photometric residuals, ARp, between 
GSC 2.3 and SDSS for a 2.5° x 60° strip along the 
celestial equator. The three panels show the resid- 
uals of 9544, 17,507, and 29,709 stellar objects 
selected within the magnitude range Rp =15.0- 
15.5, 16.5-17.0, and 17.5-18.0, respectively. The 
thick solid lines show the running means of the 
residuals. The crosses at the bottom indicates the 
position of the centers of the plates covering this 
equatorial region. 

curacy via comparisons with the GSPC-II and the 
SDSS, once the reference magnitudes have been 
transformed into the natural plate bandpasses by 
means of the synthetic color transformations dis- 
cussed in Sect. 3.3.2. 

The small-scale test shown in Figure [23] gives 
the residuals ARp with respect to SDSS for a strip 
of 2.5° X 60° along the celestial equator covered by 
a dozen plates. The thick solid lines are running 
means which evidence a zero-point systematic er- 
ror as a function of R.A., clearly correlated to the 
plate pattern. Peak-to-peak differences of about 
±0.15 mag are present, although the standard de- 
viation of the zero-point variation does not exceed 
0.05-0.06 mag. The presence of a mild magni- 
tude term is also revealed by comparing the three 
panels which show the residuals in the magnitude 
ranges Rp =15.0-15.5, 16.5-17.0, and 17.5-18.0. 



As expected, the random scatter around the mean 
clearly increases as a function of magnitude, from 
an ~ 0.10 mag to 0.14 mag of the first and last 
panels, respectively. 

An all-sky comparison is the one against GSPC- 
II, which we call semi-internal as this catalog was 
used to calibrate the plates, even though a per- 
centage of the matched objects did not partici- 
pate in the calibration because many objects were 
rejected during the pipeline reduction for various 
reasons (e.g. color outside the vahd range of the 
transformations, position close the plate border, 
etc.). Figure [24] shows the distributions of the 
residuals, ARp, between GSC 2.3 and GPSC-II 
magnitudes computed for 204,322 and 175,364 
objects down to the plate limit in the northern 
and southern hemisphere respectively. The dis- 
tribution is well behaved, with a maximum close 
to zero within a few hundredths of magnitudes 
and a standard deviation of about 0.3 mag, es- 
timated after outliers rejection at 3(7. We note 
an excess of objects in the tail of the distribu- 
tions with respect to a pure Gaussian, and also a 
slight positive skewness formed by objects having 
^Gsc 2.3 ^ ^GSPC-II ^ possibly due to unresolved 
binaries which result systematically brighter on 
the photographic plates. The residual distribution 
of the point-like objects is also shown in the right 
panel; as expected, it appears similar but with a 
smaller dispersion. 

The bottom panels show the same statistics for 
objects at galactic latitudes higher than |6| = 30°, 
where the crowding is less critical. Here the limits 
of the photometric accuracy can be directly tested; 
in the case of stellar objects, dispersions attain 
0.13 and 0.16 mag in the northern and southern 
hemisphere respectively. The results of the statis- 
tics of all the distributions shown Fig[24]are listed 
in Table [9] 

An independent analysis of GSC 2.3 was done 
with respect to the SDSS DR5 that provides 
high accuracy ugriz photometry which was trans- 
formed into the GSC-II photometric system, 
BjRpIjy, and compared with the photographic 
magnitudes. In Tables [10] , [HI and [12] we report 
the standard deviation of the residuals of stellar 
objects as a function of magnitude, in the northern 
and southern hemispheres. In the north, statistics 
are given for a 20% sample of randomly selected 
objects, while all the southern matches were used. 
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Fig. 24. — Distributions of the ARp residuals between GSPC-2 and GSC 2.3 magnitudes for all the sky 
(top panels) and high galactic latitudes (bottom panels). The solid and dashed lines indicate the statistics 
computed in the northern and southern hemisphere, respectively, for all the objects (left panels) and stellar 
objects only (right panels). 



The listed values were computed after rejecting 
3(T outliers. Figure [21] shows the same statistics 
for a subsample of 2000 objects per 0.5 mag bin, 
randomly drawn from the whole sample. 

At intermediate magnitudes, we attain a preci- 
sion of UABj — 0.16, (taRf — 0.12, and cjaik — 
0.13. The errors increase for brighter and satu- 
rated stars, as well as for the faintest object up to 
0.20-0.25 mag in the second-last magnitude bin 
before the plate limit. Note that these values, 
which are averaged on a large sky area, include 
both random and position-dependent systematic 
errors (see Fig. [23l) . 

The global zero points of the stellar objects, 
(ABj), A{Rf), A(/jv), averaged over the entire 
sky area where GSC 2.3 and SDSS DR5 overlap 
(about 8000 square degrees, mostly in the north- 
ern hemisphere) are usually very small, within a 



few hundredths of magnitude, as expected by local 
systematic errors which randomize on large scales. 
This is not true in the case of extended objects 
that, as shown in Figure [26l appear systemati- 
cally brighter in GSC 2.3 with respect to SDSS. 
The magnitude offset. Am, increases monotoni- 
cally toward brighter magnitudes up to a few mag- 
nitudes, because extended images such as galaxies 
or blends are less affected by the saturation of the 
photographic emulsion than point-like objects. In 
fact, Am -^ for fainter sources which present 
densities in the linear part of the characteristic 
function of the photographic emulsions. Figure 
1261 also shows a certain number of bright sources 
having Am « 0. Actually, these correspond to 
real point-like objects which have been misclassi- 
fied by the classification algorithm. 

Local systematic errors have been investigated 
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Fig. 25. — Stellar objects in the northern (left panel) and southern (right panel) hemisphere. Distributions 
of the Ai3j Ai?i? A/at residuals between GSC 2.3 and SDSS magnitudes for stellar objects in the northern 
hemisphere. Stars are uniformly distributed in magnitude, with a density of 2000 objects per 0.5 mag bin. 
The solid lines indicate the mean and the boundary lines at ±2cr. 






Fig. 26. — Nonstellar object in the northern (left panels) and southern hemisphere (right panels). Distribu- 
tions of the Ai3j Ai?F A/a? residuals between GSC 2.3 and SDSS magnitudes for nonstellar objects in the 
Northern hemisphere. Stars are uniformly distributed in magnitude, with a density of 2000 objects per 0.5 
mag bin. 
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Fig. 27. — GSC 2.3 - SDSS residuals. Distribution of the mean zero point, {AB.f), of stellar objects in the 
sky area 150° < a < 220° and 20° < S < 50°. The mean residuals have been computed within square regions 
of about 30' X 30'. 
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Fig. 28. — GSC 2.3 - SDSS residuals. Color map of (ARp) (left panel) and aABp (right panel) of stellar 
objects measured in the sky area 150° < a < 220° and 20° < S < 50° with a spatial resolution of about 
30' X 30'. The white crosses indicate the centers of the POSS-II plates. 
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in details within a test area of 1700 square de- 
grees selected at (150° < a < 220°, 20° < ,5 < 
50°). There, magnitude residuals were computed 
in small regions of about 30' x 30' and for dif- 
ferent magnitude bins. The left panel of Figure 
[27l shows the distribution of the mean residuals, 
(ARp), of the stellar objects in the magnitude 
range, 18.0 < Rp < 18.5, which corresponds to 
the typical magnitude limit of the GSPC-II refer- 
ence stars used for the photometric calibration of 
the plates. 

The zero-point distribution appears well be- 
haved and roughly Gaussian, although with bigger 
tails. It presents a mean of -1-0.05 mag and a stan- 
dard deviation of 0.09 mag whose value represents 
the typical scatter of the local systematic errors 
that are present in the GSC 2.3 photometry. 

The right panel of Figure [27] compares the 
zero-point distribution of three magnitude bins, 
15.0 < Rf < 15.5, 18.0 < Rf < 18.5, and 
19.5 < Rf < 20.0. The two brightest bins show 
similar distributions, having the same standard 
deviation, (T(^/^j^j^-^ ~ 0.09 mag, but a relative shift 
of a few hundredths of a magnitude. The faintest 
magnitude bin is close to the plate limit and shows 
a zero-point distribution which appears asymmet- 
ric and with a larger dispersion, (J(aRf) — 0-12 
mag, with respect to the previous ones. This de- 
pends on the lack on some plates of faint pho- 
tometric reference stars needed to fit the nonlin- 
ear density to magnitude transformation described 
in Sect. 3.3.2. Note that - 90% of the sky has 
been calibrated with photometric sequences down 
to Rf > 18 and Bj > 19.5, while for - 50% of the 
sky the photometric sequences attain Rf ^ 19.5 
and Bj > 20.5. Careful catalog users can look at 
the 4th, 5th, and 6th digits of the source status 
flag ^ in the export catalog in order to check if 
magnitudes are calibrated by means of interpolated 
or extrapolated transformations. 

Finally, the scale length and spatial properties 
of the photometric systematic errors are graphi- 
cally represented in Figure [28] which show maps 
of the photometric zero point, (ARf), and the lo- 
cal random error, a(ARF), of stellar objects with 
18.0 < Rf < 18.5 over the wide test field of about 
1700 square degrees. The patchy pattern shown 
in the left panel of Fig[28] is clearly correlated to 
the distribution of the photographic plates, whose 
centers are located on a grid having a step of 5° 



indicated by the crosses in the figure. The zero- 
point scatter is generally consistent with the 0.09 
mag standard deviation derived from the distri- 
bution shown in Fig. 1271 However, here we note 
the presence of a plate located at a « 170° and 
S « 34° severely affected by large zero-point er- 
ror because of the missing of a deep photometric 
sequence in this region. 

The right panel of Figure [55] shows the maps 
of the local random photometric error, <y/AR.F)^ 
which appears quite uniform, thanks to the selec- 
tion criterions described in Sect. 4.2.1. It is also 
interesting to note that a significantly higher er- 
ror is only present on the borders of the critical 
plate above discussed and which overlap with the 
neighbor plates. 

Similar error statistics have been found for the 
other GSC-II passbands. 

5.5. Classification 

The reliability and efficiency of the classifica- 
tion procedures in the GSC 2.3 are a complex 
interaction of many effects: plate-based imaging 
variations, external galaxy/star distributions, in- 
ternal calibrations such as object magnitude de- 
terminations, etc. A complete performance eval- 
uation is beyond the scope of this discussion and 
we limit ourselves to finding the reliability of the 
stellar classification because this is what is crucial 
for the operational uses of the GSC 2.3. 

As already mentioned, we have tested our clas- 
sification against various external catalogs; for this 
discussion we will limit ourselves to the SDSS DR5 
comparison but note that the results with other 
catalogs are similar. 

In the SDSS there are eight classifications 
which, in addition to star and galaxy, include 
cosmic rays, defects, star trails, and other arti- 
facts. We consider here only the objects labeled 
as galaxy and star. For each HTM region with 
more than ten matches we calculated the stellar 
classification success probability by comparing the 
number of SDSS stars classified as GSC 2.3 stars 
with the total number of SDSS stars and similarly 
for galaxies and nonstars. 

We limit our comparison to objects in the mag- 
nitude range Rp — 15-20 as objects brighter than 
15 are dominated by stars, and for objects fainter 
than 20 we do not have enough pixels to classify. 



35 





Fig. 29. — Left panel: ratio between the number of correctly classified GSC 2.3 stars/galaxies and the total 
number of matched SDSS stars/galaxies, as a function of Rp magnitude (completeness); Right panel: ratio 
between the number of correctly classified GSC 2.3 stars/galaxies and the total number of GSC 2.3 classified 
as stars/galaxies matched with SDSS as a function of Rp magnitude (sample purity). 
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Fig. 30. — Probability of a GSC 2.3 classification agreement with the SDSS {left: stars, right: galaxies) as a 
function of galactic latitude, and for different magnitude ranges. 
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In the left panel of Figure [29] we plot the median 
probability for the sampled regions as a function 
of F magnitude. The percentage of stars correctly 
classified drops with magnitude from better than 
90% at 18.0 to less than 50% at 20.0, while the 
probability that a galaxy is classified as nonstar in 
the GSC 2.3 remains above 90% for all the mag- 
nitude range. 

We expect a variation in galactic latitude but 
to examine this we restrict the magnitude range to 
Rp = 17 — 19 where the number of galaxies and 
stars is approximately equal and the probability 
of a correct classification is better than 80%. In 
Figure [30] we show the variation of the probability 
of a correct classification as a function of galactic 
latitude and magnitude. The probability of a cor- 
rect classification for stars drops drastically from 
better than 80% at latitudes greater than 50 de- 
grees to less than 40% for latitudes less than 20 
degrees. The galaxy /nonstar agreement remains 
better than 80% but is worse at the poles than at 
the plane. 

As previously stated, the purity of the stellar 
sample, i.e. that an object classified as a star is 
indeed a star, is important for guide stars use in 
telescope operations. This leads to the strategy 
of adopting a classifier that conservatively classi- 
fies objects as "stars" , but accurately classifies ex- 
tended objects as "non-stellar" . These properties 
are refiected in Figure [SD] the low "accuracy" of 
the stellar classifier at faint magnitudes is a result 
of the classifier being conservative, which does not 
affect operations negatively as the number of stars 
is ever increasing at fainter magnitudes. 

The right panel of figure [^H] plots the reliability 
or purity of the stellar classification as a function 
of magnitude. This is calculated by comparing the 
number of objects classified as stars both in SDSS 
and GSC 2.3 to the total number of GSC 2.3 stars. 
From this we can see that our stellar classification 
is reliable at the 90% level to fainter than 19th 
magnitude. 

5.6. Completeness 

The magnitude limit of the Schmidt plates on 
which the GSC-II is based is nominally 20.5 in Rp, 
so that the catalog is expected to be complete to 
at least Rp = 20. This completeness limit is con- 
firmed at high galactic latitudes using the SDSS 



(Drimmel et al. 2007, in preparation) along an 
equatorial strip covering more the 300 square de- 
grees; more than 99% of the SDSS point sources 
are successfully matched to a GSC 2 object down 
to magnitude Rp « 20 using a search radius of 
5" , while more than 97% of the SDSS extended 
objects were matched at magnitude Rp — 19. 
However, in very crowded fields near the Galac- 
tic plane the GSC-II may suffer incompleteness at 
magnitudes brighter than Rp = 20 due to crowd- 
ing effects. For this reason completeness in the 
GSC-II is not uniform over the sky. 

To assess the completeness limit for each HTM 
region we use a model-independent approach 
based on the assumption that, for a small mag- 
nitude interval (1-1.5 mag), the slope of the log- 
arithmic magnitude distribution, d log n{m)/ dm 
{n{m) being the object field density at magni- 
tude m, is nearly constant and positive, consistent 
with the expectation that the number of objects 
increases at ever fainter magnitudes. The con- 
stancy of the slope is clearly an approximation; 
for galaxies this assumption is very nearly true for 
sufficiently large areas of the sky, while for stars 
the slope changes slowly. In any case, these as- 
sumptions are reasonable for the magnitudes we 
are considering, whether we are considering stars 
or galaxies or all objects together. Thus, we bin 
the objects in 0.5 mag bins to magnitude Rp = 20 
for each HTM region and apply the following pro- 
cedure to estimate the completeness limit of each 
HTM region. 

First, the magnitude bin containing the maxi- 
mum number of stars is adopted as a first guess 
of the last (faintest) complete bin. This guess is 
checked by performing a linear fit to the log counts 
of the previous two bins. If this presumed last 
complete bin is the last magnitude bin, or if the 
fainter bins are empty, then a simple completeness 
test is performed: if the count in this bin agrees 
to within two times the Poisson error of the count 
predicted from the linear fit, then the maximum 
magnitude of this bin is adopted as the magni- 
tude limit. For the large majority of regions more 
than 30 degrees from the Galactic plane, the last 
bin {Rp — 19.5 to 20) meets this criterion and 
the limiting magnitude of the counts {Rp — 20) is 
adopted as the completeness limit. 

If the above criterion is not satisfied then a re- 
fined estimate of the last complete bin is made by 
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magnitude 



Fig. 31. — Example of determination of the mag- 
nitude limit. The histogram shows the raw GSC 2 
object counts from HTM region S2200131, as a 
function of Rp magnitude. The vertical arrow in- 
dicates the last complete bin, the dash-dotted line 
shows the expected magnitude distribution, the 
thick solid curve represents the actual magnitude 
distribution, and the dotted line indicates the de- 
termined 90% completeness limit. 

requiring that 1) the slope of the log-magnitude 
distribution decreases less than 30% over 1 mag, 
or that the count is not less than 2a of the count 
predicted from a linear fit to the previous two mag- 
nitude bins, and 2) that the local slope of the log- 
magnitude distribution is positive. The slope and 
change in slope of the counts are evaluated using 
a second-order polynomial that fits the counts in 
the hypothesized last complete bin and the previ- 
ous two bins. We also require that the slope of 
the above-mentioned linear fit be positive. If any 
of these conditions are not satisfied, the previous 
(brighter) magnitude bin is tested. This procedure 
is repeated until a magnitude bin is found that sat- 
isfies the above conditions. Note that if the actual 
count in the presumed last complete bin exceeds 
the predicted counts it is taken as being the last 
complete bin. 

Once the last complete bin is determined, we 
estimate the expected counts after the last com- 
plete bin. The expected counts are based on a 
local linear fit to the faintest part of the logarith- 
mic magnitude distribution determined to be com- 
plete, though now based on three magnitude bins, 
including the last complete bin if the count in this 



Fig. 32. — Fraction of sky less than 90% complete, 
as a function of Rp magnitude. 

bin does not differ more than 2cr from the esti- 
mated count based on the previous two bins. If 
the last complete bin does not satisfy this crite- 
ria then the previous three bins are used in the 
linear fit. This linear model is our expected log- 
magnitude distribution, that is log rig. 

A check is then made to see if any fainter bins 
are within 2 sigma of our expected magnitude dis- 
tribution, or exceeds the predicted counts. If so, 
the faintest such bin is considered as our faintest 
complete bin, though we retain the expected mag- 
nitude distribution already found. This check is 
necessary only to handle a few regions with par- 
ticularly "noisy" counts. Indeed, the complexity 
of the procedure described above was necessitated 
to assure that a last complete bin could be ro- 
bustly determined for all HTM regions, including 
the relatively few pathological cases that can oc- 
cur in very crowded regions where the photometric 
calibration was problematic. 

We now make a more refined estimate of the 
completeness limit using a second-order polyno- 
mial to describe the part of the log-magnitude dis- 
tribution where incompleteness sets in, from the 
bin preceding the last complete bin, and up to two 
bins beyond the last complete bin, using only bins 
with nonzero counts. Thus, since at this point in 
the algorithm the last complete bin does not coin- 
cide with the last nonempty bin, the polynomial is 
fit to three or more points. We adopt this polyno- 
mial as the measured log-magnitude distribution, 
logn. To determine to what limiting magnitude 
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Fig. 33. — Position on the sky of the HTM regions 
showing less than 90% completeness at magnitude 
Rf = 18.5. 



the counts are, say, 90% complete, we find the 
root to the equation 

\ogn{m) — logne(m) — log(.9) = 0, 

that is, the magnitude at which the ratio of the 
measured and expected counts, n/rig, is equal to 
0.9. 

Figure [3T] shows an example of the application 
of this procedure to the counts of a region showing 
evidence of incompleteness. 

Using the above procedure we estimate the 90% 
completeness limit for each HTM region using the 
total object counts. As a statistical measure of 
completeness we show in figure [32] the fraction of 
sky that is less than 90% complete as a function of 
magnitude. For example, 10% of the sky shows in- 
completeness at magnitude Rp = 18.5 or brighter, 
while at 17th magnitude 3% of the sky shows evi- 
dence of incompleteness. 

Incompleteness at magnitudes brighter than 
about Rp ~ 19 is strictly related to the object 
field density, as in this case incompleteness arises 
not from detector (plate) sensitivity, but from the 
photographic plates being unable to distinguish 
objects from the "background" light of fainter ob- 
jects whose images are blended due to the finite 
resolution of the plates. 

Figure [33] shows the location on the sky of the 
HTM regions showing incompleteness at magni- 
tudes brighter than Rp 18.5. 



6. Conclusions and future work 

In Table IIS) we compare the requirements 
provided in the GSC-II implementation plan — 
originally foreseen to reach magnitude 18 in V — 
with the actual performance of the current version 
of the catalog down to the same magnitude limit. 
In all items we have met or exceeded the specifi- 
cations except for the proper motions. 

Table [14] summarizes the main characteristics 
of GSC 2.3, which essentially superseeds the re- 
quirements in that it reaches 2-3 mag depeer. The 
astrometric and photometric errors reported in the 
table are derived from those of Tables [5] and [TTl 
and positional uncertainties are obtained by sum- 
ming in quadrature the error in each coordinate. 
A consideration of the global statistics shows that 
even with the inclusion of fainter objects we are 
able to meet the original specifications. 

Whilst no further improvements are required 
(or funded) for its use in HST operations, we 
intend to continue development to produce the 
best possible catalog for scientific as well as op- 
erational uses. Our future plans include an as- 
trometric recalibration using the UCAC 3 and 
applying a magnitude-dependent correction; this 
should not only reduce the absolute positional er- 
rors but amend the systematic error that is af- 
fecting the computed proper motions. A photo- 
metric re-reduction of all plates will also benefit 
from the improved overall quality of the available 
calibrating sequences. Finally, we plan to reclas- 
sify the objects splitting the nonstar classification 
into galaxy and blend as well as better determine 
the magnitudes of galaxies using a n algorithm that 



does not assume a stellar profile (jPetrosian et al. 
20071) . 



As of cycle 15, GSC 2.3 is the guide catalog for 
HST, and it is a reference catalog for the VLT and 
Gemini adaptive optics programs, where real guide 
stars are still preferred to laser "guide stars" ; it has 
also been provided to the Chandra, Galex, XMM- 
Newton and Swift missions. In the preparation 
for the Gaia mission, the current catalog version 
is used as a snapshot of what G aia is expected 
to observe ( Drimmel et al.l 120061) and as a base 
for the Initial Gaia Source List being compiled by 
OATo for the Gaia data reduction. Finally, GSC- 
II future release(s) will be at the heart of the As- 
trometric Support System of the very ambitious 
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Large Sky Area Multi-Object Fiber Spectroscopic 
Telescope (LAMOST) undertaken by the Chinese 
Academy of Science, and will be p art of the guide 
syste m for JWST (fspagnal 1200 ll : IStvs and Krisd 
2003h . 

In addition to the operational uses of the 
catalog itself, the GSC-II database has also 
been mined for ma ny scientific s t udies, such as 
planetary neb ulae (IKerber et al. 20041 ) . young 
open clu sters (ISciortino et alj 120001 ) . halo white 
dwarfs dCaro Uo et al. 20061) . pecuhar objects 
(iMirabelet al. 2001: ICarollo et al.l l2002l). X-ray 



pulsars (jPanzera et al.l 120031 ) . and the struc- 
ture and kinematics of stellar populatio ns of 
our G alaxy ( Vallenari et al.i i2006i : iKinman et al. 
20071) . 



Finally, it is worth going back to the den- 
sity map presented in Figure [HI as it does il- 
lustrate GSC 2.3: it provides accurate positions, 
magnitudes in three bands, and stellar classifica- 
tion for 4tt steradians to levels of completeness 
commensurate with the underlying stellar den- 
sity. Future versions extracted from the GSC-II 
II database will only improve on this situation, in- 
cluding proper motions and increasing the catalog 
accuracy as a result of improvements in reference 
catalogs and reduction procedures. 
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Table 1 
Plate material utilized for the construction of the GSC-II. 



Source^ 


Survey 


Declination 


Epoch 


Emulsion 


Band 


Depth 


Fields 


In CSC 


Code 




Range 




+Filtcr 




(mag) 




2.2/2.3 


n(i) 


Pal-QV 


(5 > 0° 


1983-85 


IIaD+W12 


l/l2 


19.5 


616 


N/Y 


S(2) 


SERC J 


.5 < -15° 


1975-87 


IIIaJ-|-GG395 


Bj 


23.0 


606 


Y/Y 


S(2) 


serc EJ 


-15° < .5 < 0° 


1979-88 


IIIa,I-|-GG395 


Bj 


23.0 


288 


Y/Y 


XE(3) 


POSS-I E 


S > -30° 


1950-58 


103aE-ficd plexiglass 


E 


20.0 


935 


N/N 


xo 


POSS-I O 


S > -30° 


1950-58 


103aO unfiltered 


O 


21.0 


935 


N/Y 


XJ 


POSS-II J 


.5 > 0° 


1987-00 


IIIaJ+GG385 


Bj 


22.5 


897 


Y/Y 


XP 


POSS-II F 


(5 > 0° 


1987-99 


IIIaF+RG610 


fli. 


20.8 


897 


Y/Y 


XI 


POSS-II N 


5 > 0° 


1989-02 


IV-N +RG9 


In 


19.5 


897 


N/Y 


xs 


aao-ses 


S < -15° 


1990-00 


IIIaF+OG590 


Rf 


22.0 


606 


Y/Y 


er 


SERC ER 


-15° < <5 < 0° 


1990-98 


IIIaF-(-OG590 


Rf 


22.0 


288 


Y/Y 


is 


SERC I 


5 < 0° 


1990-02 


IV-N -I-RG715 


In 


19.5 


731 


N/Y 


is 


MW Atlas 


(5 < 0° 


1978-85 


IV-N -I-RG715 


In 


19 


173 


N/Y 


xv(*) 


SERC-QV 


-70° < <5 < 0° 


1987-88 


IIaD+GG495 


V495 


14 


94 


N/Y 


GR*'^' 


aao-sr 


-70° < i5 < 0° 


1996-99 


IIIaF-(-OG590 


Rf 


20 


118 


Y/Y 



^Thc source code used by the GSC-II database to identify individual plates. For example, XJ442, XP442, and XI442 
identity the blue, red, and near-infrared plates of the 442th POSS-II field. 

Survey-filter abbreviations: see NOTE below. 

The Pal-QV survey was taken specifically for GSC-I. Schott filter GG495 replaced W12 starting 3 June 1984. Scanned 
with 25 micron sampling (1.7 arcsec/pixel). 

S plates originally scanned with 25 micron sampling, 440 plates at lower latitudes rescanned with 15 micron sampling 
(1.0 arcsec/pixel). 

POSS-I E filter was red plexiglass no. 2444; scanned at 25 micron sampling except 123 plates scanned at 15 micron, 
mostly south of —18°. 

SERC "Quick V" survey taken specifically for GSC-I with 4-min exposures to cover crowded southern Milky Way fields 
{\b\ < 15°, —112° < I < 34° plus two plates on the LMC). (Three similar "XX" plates were centered on M31 and each of 

the Magellanic Clouds.) 

AAO-SR survey taken specifically for GSC-II with 5-min exposures to cover crowded southern Milky Way fields 
{\b\ < 15°, -112° < / < 34° plus two plates on the LMC). 



Note. — Here, Pal-QV refers to the Palomar "Quick V" survey JLasker et al.| lHQQCD ). SERC refers to the Science & 
Engineering Research Council surveys, POSS I and II are the Palomar Observatory Sky Surveys I and II, AAO-SES and 
-SR refer to the Anglo-Australian Observatory Second Epoch Survey and Short-Red survey, and the MW Atlas is the 
SERC I/SR Atlas of the Milky Way and Magellanic Clouds (ref?). POSS-I and Pal-QV surveys utilized a 6° grid of 
telescope pointings, with minimum plate overlap, while the other surveys used a 5° degree grid. Unless stated otherwise, 
all plates scanned at 15 micron sampling (1 arcsec/pixel). Scans performed on glass copies of the S, XE and XO plates, 
rather than on the original plates. 



Table 2 
Classification parameters. 



Parameter 



Definition 



1 
2 
3 

4 
5 

6,7,8 
9-12 
13-14 
15-30 



Integrated density 

Peak density 

Semimajor axis 

Semiminor axis 

EUipticity 
2 



'^l^'^l 



■ xy 



Texture features (Malagnini 
Spike features (Lagker et aL 
Area at 16 detection thresholds 



'1983) ) 
1990)) 
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Table 3 
GSC 2.3 (July 2006). Export Binary Table Fields and Types 



Number 


Field name 


Format 


Unit 




Notes 


1 


gscID2 


Integer*4 (J) 


id 




GSC 2.3 object id 


2 


gsclID 


Character*ll (llA) 


id 




GSC 1.1 object id 


3 


hstID 


Character*ll (llA) 


id 




HST object id 


4 


RightAsc 


Real*8 (D) 


rad 




Right ascension ICRF 


5 


Declination 


Real*8 (D) 


rad 




Declination ICRF 


6 


PositionEpoch 


Real*4 (E) 


year 




Position epoch^ 


7 


raEpsilon 


Real*4 (E) 


arcsec 




Reference error on R.A.^ 


8 


decEpsilon 


Real*4 (E) 


arcsec 




Reference error on DEC^ 


14 


FpgMag 


Real*4 (E) 


mag 


Rf 


photographic magnitude 


15 


FpgMagErr 


Real*4 (E) 


mag 




Reference error on Rp^ 


16 


FpgMagCode 


Integer*2 (I) 


- 




Filter code of Rp 


17 


JpgMag 


Real*4 (E) 


mag 


Bj 


Photographic magnitude 


18 


JpgMagErr 


Real*4 (E) 


mag 




Reference error on Bj 


19 


JpgMagCode 


Integer*2 (I) 


- 




Filter code of Bj 


20 


VMag 


Real*4 (E) 


mag 


V 


[photographic] magnitude-^ 


21 


VMagErr 


Real*4 (E) 


mag 




Reference error on V^ 


22 


VMagCode 


Integer*2 (I) 


- 




Filter code of V 


23 


NpgMag 


Real*4 (E) 


mag 


In 


photographic magnitude 


24 


NpgMagErr 


Real*4 (E) 


mag 




Reference error on I^^ 


25 


NpgMagCode 


Integer*2 (I) 


- 




Filter code of In 


29 


BMag 


Real*4 (E) 


mag 




B Magnitude"* 


30 


BMagErr 


Real*4 (E) 


mag 




Reference error on B 


31 


BMagCode 


Integer*2 (I) 


- 




Filter code of B 


47 


classification 


Integer*4 (J) 


- 


Morphological classification^ 


48 


semiMajorAxis 


Real*4 (E) 


pixels 




Image semimajor axis^ 


49 


eccentricity 


Real*4 (E) 


- 




Image eccentricity^ 


50 


positionangle 


Real*4 (E) 


degrees 




Image orientation^ 


51 


sourceStatus 


Integer*4 (J) 


- 


Object processing status flag 



Piate epoch for GSC-II objects. For Tycfio 2 objects, for which Ta "^Tg, the R.A. epoch is given. 

These astrometric and photometric errors are not fomiai statistical uncertainties but a raw and 
conservative estimates to be used for teiescope operations. 

This fieid may inciude; (a) photographic V12 or Vigs from IlaD plates, (b) Vt of Tycho-2 stars, or 
(c) Johnson V from SKY2000. 

'^Bt of Tycho-2 stars or Johnson B from SKY2000 or photographic O from POSS-I 

Image classification: — "star" (i.e. point-like object) and 3 — "nonstar" (i.e. extended object). 

Morphological parameters of the same image used for the position. 
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Table 4 
gsc 2.3 bandpass codes. 



Code 


Band 


Emulsion+Filtcr 


Survey/Catalog 





Bj 


IllaJ + GG395 


SERC-J, SERC-EJ 


1 


Vi2 


IlaD + W12 


Pal-QV 


3 


B 


Johnson B 


SKY2000 


4 


V 


Johnson V 


SKY2000 


5 


Rf 


IIIaF + RG630 


ESO-R 


6 


V495 


IlaD + GG495 


Pal-QV, SERC-SV 


7 





103aO unfiltered 


POSS-I 


18 


Bj 


IllaJ + GG385 


POSSTI J 


35 


Rf 


IIIaF + RG610 


POSS-II F 


36 


Rf 


IIIaF + OG590 


SERC-ER, SERC-SR, 
AAO-R, AAO-GR 


37 


In 


IV-N + RG9 


POSS-II IR 


38 


In 


IV-N + RG715 


SERC-IR 


41 


Bx 




Tycho B 


42 


Vt 




Tycho V 






Table 5 








Source Status Flag Codes. 



Digits^ 



Flag Meaning 



1 and 2 Number of observations, some observations maybe excluded from GSC 2.3. 

3 Centroider type used: 0=barycenter, l=circular ** , 2=elliptical, 

3=FPA -I- applied barycenter ^ , 4=multicircular, 
5=multielliptical, 6=:FPA + circular, 7 FPA -|- elhptical 

4 Quality of exported Rf magnitude: O=not present, l=fit, 2=extrapolated 

5 Quality of exported Bj magnitude: O=not present, l=fit, 2=extrapolated 

6 Quality of exported V magnitude: O=not present, l=fit, 2=extrapolated 

7 Classification unanimity: 0=mixed vote, l=unanimous vote, 3=unanimous artifact 

8 Classification voters: 0=multiple 15um scans, l=One 15um scan, 

2=:multiple 25um scans, 3=0ne 25um scan 

9 Processing status: O=completed processing, 

l=object too big to be cut out on at least one plate 

10 Deblended object: 0=single object on all plates, 

l=child (deblended) object on at least one plate 

*counting from right to left; flags for Tycho and SKYMAP objects are 99999900 and 88888800 
respectively 

''All non barycentric centroiding using Gaussian fits 

''FPA = Fractional Pixel Allocation 



46 



Table 6 
gsc 2.3 global statistics 





Northern 


Southern 


All 




hemisphere 


hemisphere 


sky 


Objects 


396 700 598 


548 892 085 


945 592 683 


Point-like objects 


88 481909 


118 206 935 


206 688 844 


Extended objects 


308 218 689 


430 685150 


738 903 839 


Rp magnitudes 


354 300 890 


452 765 492 


807 066 382 


Bj magnitudes 


307 949 064 


403 822 005 


711 771 069 


Vpg magnitudes 


154 786 490 


67 947 515 


222 734 005 


In magnitudes 


245 102118 


329 725 734 


574 827 852 


Bj and Rp magnitudes 


275 035 930 


330 759 735 


605 795 665 


Bj, Rp, and Ij^ magnitudes 


210 601441 


236 097 358 


446 668 799 


GSC 1.1 objects 


9 071325 


9 598 246 


18 669 571 


Tycho 2 stars"^ 


1 201 897 


1 321 635 


2 523 552 


Sky2000/Skymapi 


348 458 


190 364 


473 943 



Without double/multiple entries (see explanation in section 4.2): 2,519,152 Tycho- 
2 stars and 175,632 SKY2000 stars 
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Table 7 
GSC 2.3 (July 2006). All-sky cumulative counts^ . 



mag 


Nobj (i?j) 


Nobj(i?F) 


Nob j (/at) 


Nobj(Bj + i?F + /w) ' 


10.0 


20 715 


35 949 


192 984 


50 794 


10.5 


34 491 


82 556 


443 489 


175119 


11.0 


54 584 


239 047 


1010 629 


564 997 


11.5 


86 600 


706 719 


2 190 201 


1519 395 


12.0 


154 228 


1 829 779 


4 415 724 


3 367 609 


12.5 


364 579 


3 994 288 


8 089 061 


6 305 430 


13.0 


1 060 240 


7 458 931 


13 612 182 


10 657 925 


13.5 


2 603 441 


12 638 946 


21619 677 


16 943 348 


14.0 


5 056 922 


20196 571 


32 990 022 


25 837 448 


14.5 


8 668 396 


31041175 


48 953 930 


38 151152 


15.0 


13 898 051 


46 366 494 


71017 843 


54 914 864 


15.5 


21361566 


67 781 378 


101286 676 


77 348 928 


16.0 


31 792 178 


97 002 324 


142 016 145 


107 002 624 


16.5 


46 052 408 


135 717 387 


196 139 738 


145 725 200 


17.0 


65119 781 


186 270 427 


267 399 004 


195 532 864 


17.5 


90 046 889 


250 652 663 


349 775 619 


254 950176 


18.0 


12 1721 348 


327112 112 


434 137 985 


319 880 384 


18.5 


161 745 847 


405 027 755 


510 834 052 


383 714 272 


19.0 


210 802 304 


489 823 020 


562 417 267 


427113 344 


19.5 


270 784 037 


583 790 690 


574 291 825 


443 187 232 


20.0 


340 876 792 


685 566 332 


574 700 787 


446 393 632 


20.5 


417 992 984 


783 784 283 


574 823 790 


446 588 640 


21.0 


497 409 326 


805 951305 


574 827 840 


446 645 408 


21.5 


573 255 946 


806 469 272 


574 827 852 


446 675 616 


22.0 


649 675 534 


806 698 629 


574 827 852 


446 689 312 


22.5 


707 695 255 


806 828 565 


574 827 852 


446 696 416 


23.0 


711708 557 


806 937 362 


574 827 852 


446 698 848 


23.5 


711722 103 


807 017 618 


574 827 852 


446 699 680 


24.0 


711736 633 


807 061276 


574 827 852 


446 699 680 


all 


711 771 069 


807 066 382 


574 827 852 


466 699 774 



Cumulative counts from magnitude to m given in the first column 

The number of objects with computed 
lated according to magnitude steps in Rf 



The number of objects with computed photometry in the 3 passbands are accumu- 
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Table 8 
GSC 2.3 - UCAC 2/SDSS: Astrometric residuals. 



stellar objects 



extended objects 



Rp 
(mag) 



N.i 



(") 



fAa 
(") 



CAS 
(") 



eA<5 



N,,. 



(^Aa 



fAa 
(") 



O'AS 



Table 9 
GPSC-2 - GSC 2.3 photometric residuals 



eA<5 



UCAC 2 






















11.5-12.0 


786 859 


0.16 


0.20 


0.17 


0.20 


81641 


0.63 


0.67 


0.59 


0.63 


12.5-13.0 


2 688 254 


0.16 


0.19 


0.16 


0.19 


438 091 


0.57 


0.62 


0.53 


0.58 


13.5-14.0 


5 467 031 


0.14 


0.19 


0.14 


0.18 


1559 410 


0.55 


0.60 


0.52 


0.57 


14.5-15.0 


9 741 299 


0.13 


0.19 


0.13 


0.17 


4 559 916 


0.48 


0.54 


0.46 


0.51 


15.5-16.0 


16 289 992 


0.12 


0.20 


0.13 


0.18 


9 553 274 


0.42 


0.48 


0.40 


0.45 


SDSS 






















14.5- 15.0 


566 319 


0.14 


0.21 


0.14 


0.16 


37169 


0.21 


0.27 


0.21 


0.24 


15.5-16.0 


878 746 


0.13 


0.21 


0.13 


0.16 


115 608 


0.22 


0.28 


0.22 


0.25 


16.5-17.0 


1194 037 


0.13 


0.22 


0.13 


0.16 


394 648 


0.18 


0.25 


0.18 


0.21 


17.5-18.0 


1 528 006 


0.14 


0.23 


0.14 


0.17 


1011007 


0.18 


0.25 


0.18 


0.21 


18.5-19.0 


1896 983 


0.16 


0.26 


0.16 


0.20 


2 759 341 


0.21 


0.28 


0.21 


0.24 


19.5-20.0 


1 608 896 


0.23 


0.31 


0.23 


0.25 


7 965 088 


0.32 


0.37 


0.31 


0.34 



Objects 



Northern Hemisphere 
Nbr (ARf) (JAR, 



Southern Hemisphere 
Nbr (ARp) aAR, 



All objects 204 322 

Stellar objects 84 996 

AU objects (high latitude) 36 714 

Stellar objects (high latutude) 25 508 



0.01 


0.28 


175 364 


0.04 


0.33 


0.03 


0.18 


60 613 


-0.02 


0.25 


0.03 


0.18 


23 977 


0.02 


0.22 


0.01 


0.13 


12 846 


-0.02 


0.16 



Table 10 
GSC 2.3 - SDSS photometric residuals ABj 



B, 



North 
n o-ABj 



South 
n (TABj 



14.5-15.0 


44333 


0.42 


23713 


0.31 


15.5-16.0 


83030 


0.19 


47614 


0.14 


16.5-17.0 


129263 


0.17 


76672 


0.14 


17.5-18.9 


167603 


0.17 


102888 


0.17 


18.5-19.0 


203192 


0.16 


114972 


0.19 


19.5-20.0 


246944 


0.18 


118327 


0.18 


20.5-21.0 


263612 


0.22 


109137 


0.22 


21.5.22.0 


150603 


0.27 


58421 


0.27 



49 



Table 11 
gsc 2.3 - sdss photometric residuals arp 



i?F 



North 

N CTAflp 



South 

N (TAR-p 



13.5-14.0 


49732 


0.16 


26838 


0.17 


14.5-15.0 


91999 


0.13 


52713 


0.12 


15.5-16.0 


146229 


0.12 


84616 


0.11 


16.5-17.0 


198843 


0.12 


117889 


0.13 


17.5-18.0 


254016 


0.13 


140616 


0.13 


18.5-19.0 


306114 


0.15 


143052 


0.15 


19.5-20.0 


222294 


0.22 


88951 


0.21 



Table 12 

GSC 2.3 - SDSS PHOTOMETRIC RESIDUALS A/n 





North 


South 


In 


N 


CTA/n 


N 


C^A/n 


13.5-14.0 


56684 


0.17 


34795 


0.12 


14.5-15.0 


120778 


0.13 


66884 


0.10 


15.5-16.0 


186497 


0.14 


106972 


0.11 


16.5-17.0 


258430 


0.15 


150240 


0.13 


17.5-18.0 


308639 


0.19 


164387 


0.21 


18.5-19.0 


143655 


0.26 


59368 


0.25 
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Table 13 
gsc ii specifications'' and performance {vnm = 18) 



Specifications 



Performance 



Reference frame 


ICRF 


ICRF 


Absolute position error 


< 0."5 


< 0."3'' 


Relative position error over a 0.°5 field 


< 0."2 


< 0."2 ^ 


Adherence to reference frame 


< 0."15 


< O.^'W 


Proper motion error (total) 


< 0.004 " year^i 


Not released 


Photometry 






Passbands 


At least 2 (1 color) 


Bj, Rp, In, Some V and O 


Magnitude error 


0.1-0.2 mag 


<0.2 


Completeness to magnitude limit 


Yes 


For 90% of the sky ^ 


Stellar classification: 


95% 


98%° 



'^source: GSC-II Implementation Plan (JLasker and Lattanzil 119941 ). Performance statistics are 
given for objects to magnitude limit V=18, as in the original implementation plan 

''From a global comparison to DR5 of the SDSS, for objects classified as stellar. The random 
contribution for extended sources is approximately 20% worse 

'^Averaged over all plates, all magnitudes 

'^Completeness is reached outside of the galactic plane; for details see Section 5.6 

"^98% of objects classified as stars are confirmed as stellar by a comparison to SDSS; see Section 
5.5 for details 
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Table 14 
GSC 2.3 Global properties 



Total objects 

Magnitude limit 

Mean epoch of positions 

Reference frame 

Astrometric reference catalogs 

Average positional accuracy '"' 

Rf < 18.5 

18.5 < Rf < 19.5 

Rf > 19.5 

Average positional precision'^''^ 

Rf < 18.5 

18.5 < Rf < 19.5 

Rf > 19.5 

Photometric reference catalogs 

Average photometric accuracy (stellar sources)*-'^; 

Rf < 18.5 

18.5 < Rf < 19.5 

Rf > 19.5 

Completeness 



945,592,683 

Bj = 22.5, Rf = 20.5, and In 

1992.5 

ICRF 

ACT + Tycho-2 

Stellar objects 

0."28 

0."30 

0."40 

Stellar objects 

0."20 

0."22 

0."32 

GSPC2 + Tycho 

0.13 mag 
0.15 mag 
0.22 mag 

>98% 



19.5 



Extended objects 

0."35 

0."37 

0."50 
Extended objects 

0."25 

0."30 

0."44 



Up to Rf=20, I > 30°('^) 



("'Accuracy referes to the combined contribution of random and systematic errors. The values are global averages 
for the SDSS DR5 sample. Note that plate-to-plate as well as north-south discrepancies can reach a few tenths of an 
arcsecond. This is mainly due to residual systematic effects, as described in Section 5.2, which will be addressed in 
the next catalog release. 

('''These values, estimated using formula [5l only show the random part contribution to the positional error. 

(°) Based on the SDSS DR5 sample. Systematic offsets of the order few hundredths of a magnitude; photometry of 
non-stellar sources suffers from systematic errors which, for very bright objects, can be as high as w 2 mag, see Fig. 



Wyoy the Galactic plane, the fraction of sky with less than 90% completeness limit is shown in Fig. 
of magnitude. 



las function 
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